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ANGULAR DISTRIBUTION OF DEUTERONS FROM 
?Be(p,d)5Be. П. 
(SUPPLEMENT) 


J. DABROWSKI AND J. SAWICKI 
Institute of Physics, Polish Academy of Science, Warsaw 
( Received June 24, 1955) 


Having finished the present paper* we have learned of the recent results of J. B. 
Reynolds and К.С. Stauding! concerning the °Be(p,d)*Be reaction for 16.5 Mev pro- 
tons. We give here the comparison of these results with our theory for ry = 5.10713 cm. 
The angular distribution is given in Fig. 5. The absolute values of the cross-section 
are given in Table II. The agreement of the Born approximation theoretical angular 
distribution with the experimental points is quite satisfactory save for the point 
O, —8?. The Butler distribution is in contradiction to experiment. However the 
discrepancy between the Born approximation theoretical and the experimental abso- 
lute values of the cross-section seems to be rather great. It seems hard to be сш 
by the neglect of the proton —?Be interaction only. 


d - 
Table II. Absolute values of A in milibarns /ster. for the ?Be (p, d) “Ве (ground state) reaction for 
5 
16.5 Mev protons. The Coulemb correction is taken into account 


Theoretical values 
for ro = O10 CMD 


Св [es 


Experimental 
valués 


0° 
8^ 16.6 + 2.8 
12° 16.8 + 2.6 


* Acta Phys. Pol. Vol. XIV (1955) Fasc. 5, 407. 
1 The authors are much indebted to dr J. B. Reynolds for letting them know his unpublished re- 


sults. 
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TEMPERATURE CHANGES AND CONDENSATION DURING 
ADIABATIC EXPANSIONS OF AIR SATURATED WITH VAPOUR* 


By Danuta SrACHÓRSKA 
Institute of Experimental Physics, Maria Curie-Sklodowska University, Lublin 
(Received August 20, 1954) 


If a dry gas of temperature Т) expands adiabatically its temperature falls, reaching 
a certain value T,. This value may be calculated with the help of well-known equations, 
provided that the ratio of initial and final volumes is known. If before expanding the gas 
contained saturated vapour at Т), condensation may occur during the expansion, and 
owing to the condensation heat we finally obtain a temperature T, which is higher than Т,. 

The aim of this work was to measure temperature changes occuring during expansion 
in a Wilson chamber. For this purpose a resistance thermometer was used made of wire 
of about 10 и diameter. This thermometer was mounted in a Wheatstone bridge system ` 
in which a cathode ray oscilloscope served as an indicator. Measurements were performed 
either in the presence of a powerful ionizing substance, or without this substance, but 
in the presence of a field which removed the ions. 

In both cases diagrams of the dependence of the temperature decrease AT on the 
ratio of final and initial volumes v,/v, were obtained. The limit of condensation is indicated 
in these curves by an abrupt change of the slope. When the ratio v/v, crosses the critical 
expansion ratio (the limit of condensation) the charactér of the oscillogram is radically 
changed. There appears a steep section corresponding to the emitted heat of conden- 
sation. The final point of this section determines the equilibrium temperature T. This 
experimental value is in good agreement with that calculated on theoretical basis. 

If we take the existence of this steep section in the oscillogram as a criterion of conden- 
sation, then from the measured values we may calculate the limit of condensation on 
ions and on neutral aggregates. The values obtained in this way for alcohol vapour, namely 
1.150 and 1.180, are in agreement with those found by other authors by means of direct 
observation of appearing drops. 

From the inclination of the mentioned section the time of the condensation could be 
estimated. For the expansion v/v; = 1.170 in the presence of an ionizing sample it is 
of the order of 0.02 sec. 


I. Theoretical introduction 


The calculation of equilibrium temperature 

If a dry gas at T, temperature and v, volume undergoes an adiabatic expansion 
till it attains volume vg then its temperature falls to a certain value Т», according: 
MED чу absurdo ek ee ыл see c шр елы ees 


* Summary of a paper published in Ann. Univ. M. Curie-Sklodowska, 7, 117 (1954). 
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to the equation 


ТЕ. ЕЕ А 

vae 3 

where k = с с, If the gas contains a vapour which is saturated at the temperature T:, 
condensation may occur during the expansion and then, owing to the heat of conden- 
sation, one finally obtains a temperature T, higher than Т. Let us assume that at 
temperature Т the remaining vapour is saturated, i. е. that condensation continues 
until it reaches a state of equilibrium. Let us further assume that the temperature 
drops first to Т, then the condensation occurs and the temperature increases to T. 
This will be justified if the expansion takes place instantaneously, ie. so quickly 
that it may be assumed that the condensation begins only after the expansion is 


ended. With these assumptions the temperature Т may be determined by the following 


method, indicated by Thompson (1898). 

Denoting the density of the saturated vapour at Т, by оу and at Т by о, and the 
ratio v/v, by x, we obtain for vapour density after expansion but before condensa- 
tion 9,/x. 

After the the condensation 


01 Q 


x La 
where Q is the heat evolved in 1 cm? of the volume of the chamber and L the heat 
of evaporation. 
Let d be the air density before the expansion and c, су and с, the specific heats 
of air, vapour and liquid respectively, 


The heat used for warming from T, to T the air, vapour, and liquid contained 
in a volume unit will accordingly amount to 


d 
= < (Т-Т), „богы Ту, (ее 
Assuming for alcohol с; = с; = 0.5 cal/gram degree, we obtain 


dc х Шъ 
Q- Sta x 99 (г. ту, 


and therefore 


gue dc + 05 о, 
x RE 


“(7--Т,) 


і 
сіп this equation is either Cor С, de on whether the expansion occures to 
a given. pressure or to a given qme 


The ratio k = CIC, for a mixture of air and vapour may be calculated from the 
Richarz equation 


пи o 1 Рі 1 i o 


о > 
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where kı, kọ are the values of this ratio for air and the given vapour respectively, 
and рі, ps their partial pressures. In this way we first calculate 7, from equation (1) 
and then T' from equation (2). 

The solution of equation (2) is based on finding such a value T, which substituted 
into the right-hand side of the equation would give for 0 a value in agreement with 
that given in tables for saturated vapour pressures at this temperature. 

The search for such a value may be simplified using a graphic method. We plot 
once for all the curve o — o (T) which gives the relation between saturated vapour 
density and temperature, and then for each expansion we plot a straight line given 
by the right-hand side of the equation as a function of T. The abscissa of the point 
of intersection gives the temperature value sought for. 

This work is concerned with direct measuring of the temperature fall — A T during 
the expansion in a Wilson chamber for various values ot the ratio 9/01, and comparing 
the observed with the calculated values of T, and T. It is a wellknown fact that in 
a Wilson chamber the formation of drops on ions begins only at a certain value of v,/v,, 
called the critical volume-expansion ratio, and if ions are absent this critical ratio 
is higher. 

It was our aim to ascertain whether the critical ratio is in any way indicated on the 
curve which represents the dependence of AT оп 4/01. 

Moreover, it was to be expected that the knowledge of temperature changes in 
time would lead to definite deductions as to the condensation process. 


II. Evidence from literature 


In literature there are to be found but few data concerning temperature changes 
in a Wilson chamber during expansion, and they are not quite consistent. 

In 1911 there appeared a paper by Chapman (1911) who measured the temperature 
at adiabatic expansion of air saturated with water (or alcohol) vapour in a 1001 
glass container. A thermocouple prepared of wires of about 24 и thickness connected 
to a galvanometer served as a measuring instrument. Аз a result of his experiments 
Chapman expressed the opinion that the temperature 7» cannot be observed, as 
the condensation heat instantly compensates the temperature depression and the 
temperature reaches immediately its equilibrium value ТІ, 

Flood (1934) published a paper on the appearance of drops of liquid in super- 
saturated vapour. In its first part the author deals with measurements of temperature. 
For these measurements he employed a resistance thermometer with a Wollaston 
wire of 8 и diameter inserted in a Wheatstone bridge, a cathode ray oscilloscope serving 
as indicator. The precision of these measurements was rather poor, as it appears from 


-- 


# 
1 The way in which Chapman calculates Ше equilibrium temperature and obtains agreement of the 
experimental and theoretical values is not clear. 
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the values given: temperature change of 50°C produced a deviation of 18 mm on the 
oscilloscope, while the thickness of the line amounted to about 1/5 of the whole 
deviation. 

Flood measured first the fall of temperature during the expansion of dry air and 
found that the expansions are approximately adiabatic up to v/v, = 1.5 (assuming 
small expansions to be adiabatic). He then measured the decrease of temperature 
accompanying the expansion of air saturated with water vapour and obtained much 
smaller values than with dry air. 

This, however, Flood does not explain as a result of condensation heat in the whole 
mass of air, but believes that the temperature depression is the same as in dry air, 
the resistance thermometer giving false indications as a result of water condensed on it. 

It is apparent from Flood's oscillograms that the temperature remained constant 
for about 1 sec. Some authors infer from this fact that for about 1 sec after the expan- 
sion the temperature in the chamber remains constant, and so does the corresponding 
supersaturation. 

Nothing particular was noticed either by Champan or by Flood in the temperature 
changes connected with crossing of the condensation limit. 

Many papers dealing with the Wilson chamber were especially concerned with 
measurements of the critical supersaturation for condensation on ions or on uncharged 
nuclei. These measurements were always based on an adiabatic expansion in a chamber 
filled with air saturated with vapour of the liquid examined. The ratio of the final 
and initial volumes v,/v, was measured, and the inside of the chamber was observed 
by means of a lens or microscope. The lowest value of v,/v, at which dew-drops 
appeared was called the critical expansion or the limit of condensation. Then the 
corresponding supersaturation, called the critical supersaturation, was calculated. 
Most of these measurements were made for water vapour. Other liquids, mostly 
organic ones, were also examined. The results of these measurements may be found 
in Tohmfor and Volmer's paper (1938). 

The values of the condensation limits for alcohol vapour are given in Table I. 


They were measured by the above method of direct observation of the appearing 
dew-drops. 


>) NES * . v 
Table I. Critical volume expansion ratio (condensation limit) В ‚ for alcohol vapour 
v; | crit. 


X-ray ionization 


positive ions 


natural ions 
present 


Authors 


negative ions 


Flood (1934) 
Volmer and Flood (1934) 
Scharrer (1939) 
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Ш. The apparatus 


Fig. 1 is a diagram of the apparatus which was employed for measuring temperatures. 
The resistance thermometer Г was made of a thin tungsten wire of approximately 10 и 
diameter (used for the production of electric bulbs). Its total resistance amounted to 
about 680 (2. The cathode ray oscilloscope used for measurements had a screen 


(0051110 - 
graph 


Fig. 1. Schematic diagram of apparatus. В — Resistance box. К — Switch. К, — Safety resistor. T — 
Resistance thermometer. P — Ring-electrode. A — Vessel with alcohol. 


of 12cm diameter, with a strong protracted luminescence, and a low frequency 
amplifier. The sensitivity, depending on the degree of amplification, could be raised 
up to 15 mm per millivolt. The sweep oscillator imparted to the light spot on the 
screen a horizontal velocity of 25 mm/sec or 50 mm/sec. 

If the Wheatstone bridge was not in equilibrium, the momentary opening of the 
key K was indicated by the appearance of a ridge on the oscilloscope (Fig 2). If the 


Fig. 2. 


bridge was brought to equilibrium by an appropriate choice of resistance R, the opening 
and closing of the key K caused no deviation of the oscilloscope line. 

The measurement was effected in the following way: after the Wheatstone bridge 
was brought to equilibrium, the key was opened and then the expansion was brought 
about. The cooling of the thermometer T'caused a sharp deviation on the oscilloscope. 
This deviation could be measured by means of a transparent millimeter scale super- 
imposed on the screen of the oscilloscope. Subsequently the key was closed. When 
the temperature reached the former value the resistance К was changed in such 
a way, as to obtain, after opening the key, a deviation equal to the one observed during 
the expansion. The change of resistance AR was of course equal to the change which 
the thermometer underwent during the expansion. 
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In this way the calibration of the oscilloscope was avoided, the degree of amplifica- 
tion could be changed, and therefore it was possible to utilize the screen of the oscillo- 
scope in a better way. During most of the measurements a temperature change of 1°C 
corresponded to a deviation of 8—10 mm on the oscilloscope. Several thermometers, 
with the coefficients 3.32 Q/deg, 4.44 Q/deg, 3.20 Q/deg, and 2.16 Q/deg were used. 
Their calibration was effected after connecting them to the bridge system, by means 
of successive immersions in water baths at 0°C and 20°C and measurement of the 
corresponding resistance changes. This was tested by measuring temperatures during 
adiabatic expansions of dry air in a 30 liter glass flask. The differences of calculated 
and measured temperatures were of the order of 0.1°C. The last thermometer 
(2.16 Q/deg) was calibrated only with the help of expansions of dry air. 

Fig. l represents a diagramatic cross-section of the chamber used for the measure- 
ments. The diameter amounted to 15 cm, the depth to 6—7 cm depending on the 
initial position of the piston. This position determined the magnitude of the expansion. 

The expansion was effected by opening the valve connecting the space behind 
the piston to a vacaum container where the absolute pressure was of the order of 
10 mm Hg. The piston then moved violently back until it hit the wall of the chamber. 
This was therefore a volume-defined expansion. The time of the piston movement 
was of the order of 0.01 sec. 

The magnitude of the expansion was determined for every measurement in the 
following way: if before the expansion the pressure in the chamber was equal to the 
atmospheric pressure b, and after the expansion and equalization of temperature 
to b —p, and if the vapour pressure at the given temperature was x, then the expansion 
ratio amounted to 


^» b-m 
v D р я 


The air which was introduced into the chamber first passed through a CaCl, water 
vapour trap and a cotton filter of 30 cm thickness. The measurements were performed 
alternatively in the presence of an electric field of 60 V/cm, or in the absence of 
such a field, but in the presence of an ionizing sample (15 mg of radium contained 
in platinum needles placed at a distance of several cm from the chamber). 

The initial temperature T, was measured with a precision of 0.1°C by means of 


a thermometer placed in the chamber. During all measurements the initial tempera- 
ture was between 19°C and 20°C. 


IV. The measurements 


Two fundamental types of oscillograms тегей рип; They are diagramatically 
represented in Figs. 3 and 4. The horizontal line represents the path of spot on the ` 


screen of the oscilloscope when the temperature of the thermometer is equal to the 
initial temperature T}. 
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At the moment indicated by point А the expansion occurs accompanied by sudden 
drop of temperature, which causes the light spot to rise up to point B. Therefore the 
difference in height of points 4 and B measures the maximum temperature drop regis- 
tered by the thermometer. 


B 
А B 
D 
С 
ud À А 
Fig. 3. Fig. 4. 


In the oscillograms of the first type (Fig. 3) a gradual rise of temperature begins 
at point B, until the temperature of the surroundings is reached. This is indicated 


Fig. 5. Fig. 6. 
Fig. 5. Oscillogram of an expansion v/v, = 1.168 with sweeping field on. The photograph was made 
through a transparent plate with a mm scale which is easily seen in the upper part of the curve 
| The segment at bottom left represents 1 sec. 

Fig. 6. Oscillogram of an expansion v,/v, = 1.170 taken in the presence of an jonizing agent. Sweeping 
_ field off. The photograph was made trough a transparent plate with a mm scale which is easily seen in 
the upper part of the curve. The segment at bottom left represents 1 sec | 
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by the gentle slope of the curve BC. Figs 5 and 8'are photographs of oscillograms of 
this type. 

On oscillograms of the second type (Fig 4) one first observes, as before, a sudden 
drop. of temperature; following that however, one observes a rapid rise up to the 
level of point D; only then begins a gradual rise up to the surrounding temperature. 
Figs 6, 7, and 9 are photographs of this 
second type of oscillograms. 

There also occurred oscillograms of an 
intermediate type, i. e. such in which no 
breaking point occurred, but the rapid rise 
of temperature passed continuously into 
a slow rise. 

The occurrance of these types of oscillo- 
grams depends on the magnitude of the ex- 
pansion. For small expansions one obtains 
the first type. When they gradually become 
larger, the value is reached at which the inter- 
mediate type appears. When the ratio v/v, 
is slightly higher one obtains an oscillogram 
distinctly of the second type. 

At stil larger expansion the breaking 
point becomes more and more distinct. The 
steep section gradually approaches a vertical 
position. It also becomes gradually shorter 


for two reasons: in the first place, because 
the breaking point takes up higher positions 


Fig. 7. Oscillogram of an expansion v,/v, = 
1.178 taken in the presence of an ionizing 2 
agent. Sweeping field off. The segment ай and in the second place because the total 


bottom left represents 1 sec. height of the oscillogram does not rise in 


proportion to the magnitude of the expansion, 
but either rises very liitle, or even within a certain range diminishes. 
In some cases, to which we shall return later, the vertical section disappears totally 
and we obtain, with large expansions, an oscillogram of the first type. 

The magnitude of the expansion at which the first type turns into the second 
one is smaller in the presence of ions than when we remove the ionizing substance 
and apply an electric field. Oscillograms on Figs 8 and 9 were obtained for the same 
expansion ratio v/v, = 1.170. In the first case, however, the ions were removed by 
means of an electric field, while in the second case (Fig. 9) the field was absent and 
the ionizing sample was placed near the chamber. It is obvious that at this expansion 
the presence of a large number of ions causes the appearance of the steep section. 

Figures 10, 11 and 12 present the results of three series of measurements with 


vapour of absolute alcohol. The conditions were the same for all three series. Similarity 
of the results proves their reproducibility. 


va 
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The thick broken line is а diagram of tem 


perature changes in the presence of the 
electric field which removed the ions. 


The thick continuous line is the corresponding 
diagram in the case of strong ionization. 


If at a given measurement the oscillogram belonged to the first type, then on the 
diagram this measurement is indicated by one point, the coordinate of which corresponds 


Fig- 8 Fig. 9 


Fig. 8. Oscillogram of an expansion v/v, = 1.170 with sweeping field on. The curve shows small oscil- 
lations of frequency 50 cps. 
Fig. 9. Oscillogram of an expansion v,/v, = 1.170 taken in the presence of ап ionizing agent. Sweeping 
_ дей off. The curve shows small oscillations. of frequency 50 cps. 


to the maximum temperature drop. If the oscillogram was of the second type, then 
on the diagram we have two points connected by a broken. vertical segment. The 
upper point corrensponds to the maximum point B (see Fig 4), and its ordinate 
indicates the maximum temperature drop registered. The er point corresponds 
to the point D, and its ordinate represents the end of the rapid rise of erei 
after expansion. The intermediate type is represented by one point at which a broken 
curve starts downwards, without, however, another point to terminate it. 

| Besides these experimental curves we plotted on the same diagram two curves 
obtained from theoretical calculations. The uppermost curve p presents the МЫ ГИ 
| change calculated from equation (1), ie. the change which would take place по 
| ` condensation heat were evolved. The curve s gives the changes оса from nue 
{ (2), i.e. the changes which would take place if complete condensation occurred up 
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to the point of equilibrium between liquid and vapour phases. It was assumed in 
this case that the heating of air due to condensation takes place at constant volume, 
that is in equation (2) the specific heat of air is equal to c,. 

It is striking that temperatures corresponding to breaking points of the oscillo- 
grams should lie on curve s. This may be explained if we assume that the steep section 


AR 
—о— ionizing agent present, field off 
AT 
à А +„—@+— 407209 agent absent, field on 
.—- oscillograms of the transient type 
1 
-Ф- oscitlograms of the second type 
20 
15 = 


6 
| Ё 
| : ت‎ 
Fig. CREUSE, ues of temperature drop ж à шев ofthe expansion ain’) =: 
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expansion corresponding to an oscillozram of the first type. With the next measuring 
At is obvious that the limit of condensation must be placed between the values of 
the abscissse of both these points? 
Denoting by E, the limit of condensation on ions and by E, the limit on uncharged 
nude we obtain 


— == сбит орет present „Пеша of 
se UN обет Gsen, field ол 
Om sp С50йорпт of the transient fgpe 


““О-- е-2»-- 0520005 of the зегоаа {ge 
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On comparing these values with those given in the table Гоп р. 8 it is seen that 
the value of Е; found by ше is lower than those quoted concerning condensation on 
natural ions. This is of course explained by-the strong ionisation caused by the radium 
sample. 

In the literature on this subject (Tohmfor and Volmer 1938, Frenkel 1945) experi- 
mental and theoretical data are given, proving that the condensation limit shifts 
in the direction of lower expansions as the number of ions increases. 


AR —о— тилу agent present, field off 
p 
p |-47 „—@»— (0Пігі70 agent absent, field on 
—Q— %-4»- ascillograms of [he transient type 
РА سن‎ <>- oscullograms of the second type 


Fig. 12. Measured values of temperature drop as a function of the expansion ratio. 


With artificial ionization by X-rays Scharrer obtained for positive ions the value 
1.144 and for negative ones 1.160. The value obtained by me lies between these two. 
An exact comparison is impossible, as no data are given as to the intensity of 
ionization. 

The condensation limit on uncharged nuclei E, in my measurements had a higher 
value than those in Table I, but within the limits of error it is in agreement with the 
value 1.178 as given by Flood. 

Changes in the oscillogram with the growth of expansions indicates the way in 
which the condensation velocity changes. This is easily seen from the inclination 
of the steep section of oscillograms of the second type. In the closest vicinity of 
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the condensation limit we have the transient type, which corresponds to a relatively 
low condensation velocity. Shortly after the limit is crossed, however, the section 
becomes very steep, so that the condensation takes place during about 0.02 sec. 
With a further increase of expansion this time becomes still shorter, but the oscilloscope 
employed by me did not permit of a more exact determination. 

The horizontal section plotted at the lower part of the oscillograms reproduced 
corresponds to l sec. ў 

On the oscillograms іп Figs. 8 and 9 there are visible current vibrations (50 cycles 
per sec), which facilitate the estimation of time. 


V. The influence of impurities 


Fig. 13 represents a diagram of several series of measurements performed during 6 
weeks without opening or cleaning the chamber. Alcohol remained in the chamber 
all the time and under its prolonged action the rubber seals partly dissolved. At the 
bottom of the chamber under the cup there appeared a drop of such a solution. 


AR 
2 


20 


15 


10 


V2 
У, 


1,05 1,10 1,15 120 1,25 


Fig. 13. Temperature drop in the presence of contaminations. Measurements relating to oscillograms 
of transient type are shown by arrows. 


Curve (1) gives the results of the earliest series of measurements, when the effect 
of alcohol on rubber was not yet visible. Next came curve (2) and then curve. (3). 
‘After the chamber was ventilated for several hours with a stream of air the broken 


curve (4) was obtained. 
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From this it is evident that impurities may cause a depression of the curve, 1.е. 
an increase of the amount of condensed steam, and besides shift the limit of condensa- 
tion in the direction of lower expansions. 

Parts of curves (2), (3), and (4) corresponding to large expansions are situated on 
the equilibrium line s. In the corresponding oscillograms the vertical section disap- 
peared, which means that for large expansions we obtain again oscillograms of the first 
type. This may be explained by the fact that now condensation takes place with great 
velocity simultaneously with the expansion and therefore a point placed on the 
equilibrium line is directly reached. 

Small amounts of a body soluble in the given liquid, when dispersed in air, may 
become condensation nuclei. This is a well known phenomenon, especially important 
in meteorology. The influence of impurities on the shape of the curves may be com- 
pletely explained by the presence of these additional condensation nuclei. 

In the series of measurements, the results of which were discussed in the former 
chapter (Figs. 10, 11, 12) special measures were taken to prevent the presence of 
the impurities. In the first place all parts of the chamber were carefully cleaned, 
then the chamber was assembled and closed. Finally the chamber was ventilated 
during over 10 hours with a stream of dry air filtered through a cotton wad. 
Following that, alcohol was introduced into cup 4 by means of a pipette, and 
after the vapour became saturated measurements were made. These measurements 
lasted at most two days. After such a series the chamber was again ventilated until 
the alcohol completely evaporated. The ventilation was then still continued for several 
hours. Then alcohol was again introduced into cup A, and the next series of measure- 
ments began. 

Curves obtained in this way are almost identical. Small differences may be ascribed 
to those impurities, which remained in spite of using the measures described. 


VI. Adiabaticity of expansions 


Measurements were made with dry air with the aim of testing to what extent the 
expansions are adiabatic. Tables II, Ш and IV give the results. The deviations from 
adiabaticity in the case of a 30 1 flask cause temperature differences of the order 
of 0.1°C. In a chamber of 500 cm? volume the difference is about 1.5?C if the expansion 
is not too large. For the volume of 1000 cm3 the temperature measured is higher on 
an average by 0.5°C than the calculated value. This is due to thermal exchange with 
the surroundings. This exchange takes place more and more slowly with the 
decrease of the ratio surface/volume,: that is with the increase of volume of the 
container. | 

These are the results with dry air. If alcohol is present in the chamber the case 
_ is somewhat different, as the breaking points of the oscillograms are then situated 
on the equilibrium line and not half a degree below it. This is due to the fact that if 
we have a suspension of liquid drops in equilibrium with the vapour then most of 
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the heat from outside is used to evaporate the droplets. The increase of temperature 
is in this case so small, that it does not surpass the experimental error.? 


Table II Dry air expansion in a 30 1 flask 


P. AT calculated А Т observed к Difference 
2 

1.102 7.9°C 7.8°C + 0.1°С 
1.105 8.2 8.3 | От 
1.061 4.97 4.90 +0.07 
1.086 6.84 6.82 +0.02 
1.075 5.93 5.84 +0.09 
1.094 7.36 7.43 —0.07 
1.120 9.47 9.37 +0.10 
1.139 10.75 10.38 -- 0.37 
1.047 3.81 3.69 + 0.12 


3 Let us consider І cm? of the chamber volume and let us assume that in it we have a spray of alcohol 
drops in air and alcohol vapour. The system is in equilibrium at temperature tọ. Let us supply the system 
with a small amount of heat — q. Аз a result the temperature will reach tọ + At, a certain amount Am 
of the liquid will evaporate, and a new equilibrium will result. 

If o(t) signifies the\saturated vapour density at temperature 1, then 


di 
Am = 0( + At) — o(t) = (4) -At 
= 


We therefore obtain 


do de | 
а= M: At + Am І = М: А: + НЕ Е B B 


і-% dt r=, 


where M signifies the thermal capacity of 1 cm? of air together with the vapour and liquid contained in it, 
and L is the heat of vaporization. 


Therefore ME q 


de 
SEL 
фа dt 
Let 4 signify the amount of heat, necessary to heat 1 cm? of air by 0,5°C, or 4 = 0.5 M. Taking 


do 


L = 930 cal/g, М = 0.0013 X 0.2 cal/deg, À L 5 X 107tgjdeg 


we obtain At = 0.1°С 
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Table IV Dry air expansion in a Wilson chamber of 1000 cm? volume 


з AT calculated AT observed Difference 
0 jt 
1.041 47°С 4.1°С + 0.6°С 
1.049 5.6 5.0 -4- 0.6 
1.042 4.9 4.5 + 0.4 
1.077 8.6 8.2 —+ 0.4 
1.077 8.6 8.0 -- 0.6 


ГП. Temperature changes at low expansions 


The part of the curves corresponding to expansions which are lower than the critical 
expansion lies above the equilibrium line, but well below the line relating to such 
temperature changes as would take place if no condensation occurred. АП the three 
curves (Figs. 10, 11, 12) representing the results of measurements in the chamber 
are in complete agreement in these parts. They are also in agreement with values 
found for alcohol by Chapman (1911). 

The fact that in this case the measurement gives a temperature which is much 
higher than that calculated from (1), may be explained if we assume, according 
to Flood, that the thermometer becomes a condensation center as soon as its tempera- 
ture falls below the dew-point, in spite of the fact that condensation does not take 
place in the whole mass of the gas. The heat evolved during the condensation heats 
the thermometer. The amount of heat which the wire receives in a unit of time is 
determined by the velocity of diffusion of the vapour. At the same time the thermo- 
meter loses heat to its surroundings. Equilibrium is reached at a certain temperature 
which is recorded by the observer. 


The following table of temperatures for the expansion of 1.140 may serve as an 
example. 


Initial temperature 20°C 
Dew-point after expansion 17.7 
Temperature calculated from eq. (1) 6.58 
Temperature observed 14.5 


The temperature indicated by the thermometer is determined by the diffusion 
coefficient of the vapour and by the thermal conductivity of air. This explains the 
agreement of results obtained by various research-workers in different conditions. 
If condensation takes place in the whole mass of the gas, then the diffusion of vapour 
towards the wire ceases because concentration differences disappear, and the thermo- 
meter soon acquires the temperature of its surroundings. To confirm this explanation 
measurements were performed using the same thermometer heated to a temperature 
higher than that of its surroundings. For this purpose the current flowing in the 
- circuit was increased to such an extent that in spite of the temperature depression 
due to the expansion, the temperature of the thermometer remained somewhat higher 


ние. m 
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than the dew-point. The points obtained from these measurements are in fact placed on 
the line calculated directly from equation (1). (Fig. 11, points indicated by asterisks) 
This seems to prove the soundness of the above explanation. 

À certain objection could be made, The thermometer in this case has a temperature 
higher than that of its surroundings by approximately 10°C. In the nearest vicinity of 
the wire, therefore, the vapour is unsaturated, while at larger distances a state of 
saturation exists. Can we infer from the fact that in the direct vicinity of the wire 
no condensation heat is evolved, that the same is true in the whole mass of the gas? 
The answer is suggested by the shape of the oscillogram. If а large temperature drop 
appeared only in the nearest vicinity of the thermometer, and in the whole mass 
of the gas there existed a temperature higher by several degrees, then after expansion 


. 


the wire would warm up more quickly and therefore the oscillogram in the presence 
of alcohol would show a more rapid decrease than in the case of dry air. 

Such a difference, however, was not observed. The inclination of the oscillogram 
is the same for expansion of dry air and of air saturated with alcohol vapour, provided 
the temperature depression is equal. Hence we may conclude that in a Wilson chamber 
during expansions under the condensation limit the temperature falls to a value which 
is near to the value calculated from equation (1). On the other hand above the conden- 
sation limit the equilibrium temperature quickly sets in. 


ИІП. The results 


To sum up, it may be stated that there are no sufficient grounds to consider the 
results of measurements for expansions under the condensation limit as true tempera- 
ture of air in the chamber, in spite of the good agreement of the results obtained 
by Chapman, Flood and the present author. 

If the condensation limit is crossed, however, it may be stated on the basis of an 
analysis of the oscillograms that the condensation does actually take place and that 


-it continues until a state of equilibrium between the liquid and vapour is reached. 


The method employed facilitates the measurement of equilibrium temperature with 
a relatively high precision. | 

It was also possible to determine the critical volume-expansion ratio for ions and 
neutral aggregates. It is interesting that when employing a criterion of condensation 
based on measurements of temperatures the results obtained are in agreement with 
results of measurements based on a method of direct observation of appearing 


` dewdrops, though the results in the first case are not, as a matter of fact, more 


exact than in the second one. M VALE T, à 
When condensation occurs it is possible to estimate the time period of its duration. 


This period becomes shorter with growing expansions. During an expansion v/v, = 
— 1.170 in the presence of an ionizing substance it amounted to about 0.02 sec. 
No difference of time was observed between condensation on ions and on neutral 


aggregates. 
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The end of the condensation and a rise of temperature of course remove the state of 
supersaturation. "Therefore in the conditions described supersaturation lasted for 
about 0.02 sec. 

Only in the case of expansions corresponding to an oscillogram of the transient 
type, that is just above the condensation limit, supersaturation could exist for a longer 
time. 

For comparison's sake we present photographs taken from Flood's paper which 
was quoted previously (Fig. 14). These photographs shows that after the expansion 


v 
Fig. 14. Dry air. ® = 116 Dry air, 2 = 138 Air with water vapour, bn 1.45 
р 


1 
] cm corresponds aproximately to 1 sec. 


of dry air, temperature remained constant for about lsec. This may be explained 
by the fact that Flood's apparatus was made of glass only, and therefore thermal 
exchange with the surroundings took more time than in a chamber possessing many 
metal parts. 

The last oscillogram gives temperature changes during the expansion of 1.45 
in air saturated with water vapour. This expansion is considerably larger than the 
critical volume-expansion ratio, which for water vapour amounts to 1.275, according 
to the quoted paper by Flood. Condensation must have taken place therefore not 
only on the thermometer, but in the whole volume of the chamber. On comparison 
with photographs presented in Figs 6, 7, and 9 it seems that the horizontal section 
in Flood's oscillogram should be treated as corresponding to equilibrium temperature 
established after the condensation has ended. The ridge I have marked with an arrow 
might correspond to the steep section on the oscillograms in Figs 6, 7, and 9. It 
could not appear more distincly as Floods electric system was not suited to register 
very rapid changes.* | 

In view of the above the conclusion from Flood’s photographs that supersatura- 
tion exists in the chamber for about 1 sec seems unjustified. 

It is a pleasure to thank Professor Stanislaw Ziemecki, Head of the Institute of 
Experimental Physics, Maria Curie-Sklodowska University, for his friendly leadership, 
the interest he took in the development of this work, and the creation of the best 
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I thank once more Professor A. Goldschmied, Dean of the Medical Faculty, Medical 
Academy, Lublin, for the loan of an electrocardiograph which facilitated the completion 
of the described experiments. To Tadeusz Lewandowski, Head of the Mechanical shop 
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* A capacity of 150 uF connected in parallel to the entrance of the oscilloscope. 
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КРАТКОЕ СОДЕРЖАНИЕ 


Стахурская, изменения температуры и конденсация при адиабатических 
расширениях воздуха, насышенноһю паром 


Если сухой газ при температуре T, адиабатически расширится, то темпера- 
тура его падет до некоторого значения Т», которое можем определить при по- 
средстве известных уравнений, зная взаимоотношение первоначального и оконча- 
тельного объёма. 

Если до расширения газ содержал пар, насыщенный при температуре Т, то 
расигирение может сопровождать конденсация M, вследствие выделения конден- 
сационной теплоты, получим в конце концов некоторую температуру Т;, выше 
чем Т. 

Задачей работы было измерение перемен температуры, возникающих в ка- 
мере Вильсона во время расширения газа. Для этого употреблено резистивный 
термометр, изготовленный из проволоки толщиною около 10 ц, включенной 
в мостик, в котором индикатором был катодный осциллограф. Измерения были 
проведены либо при сильно ионизирующем препарате, либо без препарата, но 
в присутствии ионоочистительного поля. 

В обоих случаях получены диаграммы зависимости падения температуры 
ДТ от отношения окончательного объёма к первоначальному 2 . Конденса- 

1 
ционный порог отмечается на этих кривых изменением наклонения. Если отно- 


шение da переступит порог конденсации, TO изменится в основе и характер 
v. 

бак Та. Является именно в ней крутой отрезок, соответствующий выде- 

лению конденсационной теплоты. Конец этого отрезка определяет температуру 

равновесия T, причем получаем хорошую согласованность с величиной вычислен- 

ной теоретически. 

Если появления в осциллограмме крутого отрезка считать критерием кон- 
денсации, то по произведённым измерениям можно определить порог конденса- 
ции на ионах и порог конденсации на нейтральных агрегатах. Найденные таким 
образом для алкогольного пара цифры 1,150 и 1,180 сходятся со значениями 
встречаемыми в литературе, полученными методом непосредственного наблюде- 
ния за обоазованием капель. 

По наклонению вышеуказанного крутого отрезка можно было определить вре- 


v 
EE 
мя продолжительности конденсации. Для расширения — = 1,170 в присутствии 
P 


ионизирующего препарата время это будет ряда 0,02 сек. 
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Exact formulae are given for an electromagnetic field radiated by a dipole situated 
inside а wave guide of an arbitrary cross section, the moment of the dipole depending 
on time through a factor 1(t) (1 — ее) sin wot, where 1(t) is Heaviside's step function. 


1. Introduction 


In this paper we investigate the following problem: a dipole with a moment M 
is placed inside a wave guide of an arbitrary cross section. M depends on time through 
a factor l(t); p(t), 1(t) being Heaviside's step function. We are interested in the electric 
and magnetic field intensities at an arbitrary point inside the guide, and at an arbitrary 
moment £ > 0. The cross section of the tube is supposed to be constant. 

The solution of the problem,.as formulated above, might be of some interest 
from the point of view of some applications of wave guides in radar, microwave trans- 
mission, etc., namely in all cases where the time constant of the wave guide is of the 
same order of magnitude as is the pulse duration. 

The transients in wave guides have already been investigated by a certain number 
of authors. 

‚ Соне (1945) and Rubinowicz (1954) investigated a particular solution correspond- 
ing to a field whose axial component behaved in a certain cross section as 
1(0) cos (wot — о). The paper by Rubinowicz contains a full discussion of the results 
which are presented in the form of convergent series and asymptotic expressions. 
These results provide a basis for the present paper. Other authors as Cohn (1953) 
and Gajewski (1955) evaluated numerically the appropriate functions, thus obtaining 
a physical picture of the transient process; Poincelot (1954) succeeded in finding the 
time-constant of a wave guide. The problem of a transient field due to a dipole with 
a given current was formulated by Kahan and Colombo (1948) who proposed a method 
of evaluating the scalar and vector potentials but did not present any effective solu- 


tion. 
(25) 
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2. The field at the wave front 


We consider a doubly infinite wave guide, its walls being perfectly conducting. 
The constants of the medium inside the guide are £ =u = 1, с = 0.1 

We choose a rectangular coordinate system with the z-axis parallel to the walls. 
The moment of the dipole is of the form 


M = M: f(t) 
where f(t) is the function through which the moment of the dipole depends on time. 
This function contains Heaviside’s step function and may be represented as 


f (t) = 10) -р(0). 


The current density corresponding to such a dipole can be expressed as 


dM Я 
M ile) = м (n) fO, (2.1) 
where ó(r) is Dirac's 6-function (in this section we suppose Ше dipole to be situated 
in the centre of coordinates). 
We note first (Rubinowicz — in press) that the initial disturbance propagates 
along the guide in exactly the same manner as in the case of a dipole situated in free 


j (r,t) = 


Fig. la Fig. 1b 


space. In other words, the field intensities on a spherical surface À of radius ct drawn 
inside the guide are the same as those on a sphere B, of the same radius, drawn around 
a dipole placed in free space (Fig. 1). But these last intensities are easy to evaluate 


from the known formulae 


E = — grad У — т (2.2) 


Н = rot А, (2.3) 


where the vector potential А is given аз 


T. (2.4) 


! Gaussian units are used throughout. 
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Inserting (2.1) into (2.4) we find 4. The scalar potential V can be found from the condi- 
tion that 

1 28V 

— — + div A = 0. (2.5) 


в 8: 


Substituting A and V into (2.2) and (2.3) we obtain the following formulae for E and H 
at a distance К from the dipole 
Ag 
[EE л 


и: 


Е = gg В(и-Мим-м CRE [3 (M ` N) N — M] 
za 
Mas DE и, 20 
, R 11 R 
i [rts] (8) (2.7) 
HE. ee. VE quom MXN, 
where 
wat, poc. f" (9-4. 


We now find the conditions which must be satisfied by fin order to ensure the unique- 
ness of the solution. 

It was shown (cf. Rubinowicz 1955), that an electromagnetic field can be disconti- 
nuous on a surface moving with the velocity of light, but it is necessary that 


АН Хп -Vz ДЕ, 
(2.8) 
n ^+ Ан, 
и 


ХЕ = Е, — Er ДН = H, — Hi. 


where 


E,, H, and Е», Н, are the field intensities respectively on the front side and rear side 
of the moving surface; п denotes a unit vector normal to the surface and pointing in 
the direction of the motion. 

In our case the moving surface is simply the front of the first disturbance, thus 
n — N and E, = Н, = 0. ДЕ and ДН are therefore equal to the functions (2.6) 


2. R 
and (2.7) taken at the wave front, ie. with inserted t = + Hence 
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Пази N) N — M] 
+EO Bm FO мх (N X М), (2.9) 


ше 


Substituting (2.9) into (2.8) with e =u = 1, we get 
Оо (2.10) 


These are conditions to be satisfied by the function f. 
It may be worth while to note that e.g. the functions 


f(t) = 10) sin о), ` f(t) = Ц) cos wot (2.11) 
do not satisfy (2.10). 
An example of a function f satisfying (2.10) is 
f(t) = Ц9 (1 — e™) sin ор, (2.12) 
because in that case | 
f'(t) = Цдре“ sin cxt + cy (1 — е-“) cos wat], (2.13) 


and thus the conditions (2.10) are satisfied. 

In view of these-results we shall examine the case where the moment of a dipole 
placed inside a wave guide depends on time through a factor 19 (1 — е“) sin wot. 
The corresponding current density is, according to (2.1), 


d 
‚= FM д = mr le) 10 [gee sin yt + wo (1-26) cos ox 


= Лазо (r) 1(t) Е еге sin Wot + (1 — e tt) cos он |. (2.14) 


where ds, is the *length* of the dipole, and J is the asymptotic current amplitude. 
Because of the linearity of the field equations, the resulting field due to a dipole 
with such a current density can be expressed as 


1 2 o 3 
Eps = Е — Е + >E, 
% 
i 2 1 (2.15) 
Н,„ = H — H + ÊH, 
Wo 


UTP натыз 
where E, H, Е, H, Е, H are field intensities corresponding to the following current 
densities 
= J 446 (г) 1(ї) cos wot, 
ول‎ = J ds, б(ғ) 1(t) e~ cos оу), (2.16) 
= J ds, (r) 1(t) e-*' sin wot. 
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We shall use eomplex representation and investigate a field corresponding to the 
following current density 


J = Лазо ó(r) 1 (t) ei @ i-o, (2.11) 


(we consider as physically significant the real part of this and similar quantities). 
We have put here 


^ 


Let E (00, 0,4) and И (оо, 0,2) designate the field intensities of the respective 
field. We then have 


1 Й 1 3 

Е = E (09.0.0), H === Н (0,0,0), 

2 ? 2 h 

E = Е (w,,0,)0, Н = Н (0.0.0), (2,19) 
3 y 3 + Ê 

E = E (09,012), H = Н (w,,0,7/2) 


Thus the problem is now reduced to that of evaluating the fields Ё and À. Knowing 


ita 292-623 
them we obtain E, Н, Е, Н, E, Н from (2.19). Substituting these vectors into (2.15) 


we find the field intensities of the resulting field. - 


3. Solution of the problem 


Any electromagnetic field inside a wave guide can be expressed in terms of fields 


corresponding to various modes: 
ЕЕ’. 
НЕ» He. 


Тһе vectors E and H™ can be resolved into transverse and longitudinal compo- 


(3.1) 


nents, so that 


ве gp 1n E. ЖЕ 
uo? — ну” оп H9, : 


where n is a unit vector directed along the z-axis. In a region without currents the 
transverse components Eo and H ni can be expressed in terms of the longitudinal 


components as 


ЗЕ (т) d (m) 
E m E3 grad, 22 Б cK. 2 n X grad, 2 , 
2 H” 1 à Eo 
H та gradg "SEA p pn X grads жега з (3.3) 
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9 д З 
where с is the velocity of light in vacuum, grade = тт + J 2у ‚ Ки and K, are eigen- 
values (arranged in ascending order of magnitude) of the following two problems: 
DEAL T uw | | 
те да + ko = 0 inside О, (3.4) 
ф=0. on С; 
and 
2 2 
E Е e Ф LK*$—0 inside Q, | 
NN (8.5) 
= оп С 
2v — i 


Fig. 2 


where О denotes the cross section of the tube, C is the boundary of Q, and v is the 4 
outward normal to С (Fig. 2). We note the well known fact that 


ES 
~ 
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с 


where 


Ф„ — normalized eigenfunction of the problem (3.4) belonging to the eigenvalue k 
(in case of degeneracy several eigenvalues may be equal), 
Ф„ — normalized eigenfunction of the problem (3.5) belonging to the eigenvalue A 
Ри, = 27M · grad; Pm (хо, уу), 
On = 2r[M X grad; Ф, (xo, ¥o)] ` п, (3.9) 
Ry, = 2nM «n ф„ (ху Yo). 


and 
un RE exp m i [ot — оо, || 
(242 2 >. 
Vo? — о | 
— ic 
Um) — ex а: wt — مل‎ 0? |\, 
yo? — 0,2 Е | | с ч 


the lower suffix after u” and 17°? denoting respective differentiation.? 
According to (3.3) we thus have for the steady state 


E = сы (Ри иш” — Ба? Вы из”) grad, Фи 


1 sm 
Eb ai Uf” п x grad, Dm 
1 sim 
Hy = — тту; Om Па? grad, Py (3.10) 
1 


(PLU — King? Ка uf?) n x grad, фи. 


EN XE 


The formulae (3.8) and (3.10) fully determine the steady-state field inside the guide. 
We suppose now that the current in the dipole depends on time not through 


-il t—a) 
> 


the factor е as previously assumed, but through the factor l(t) e 
where 0% is, in general, a complex number. Тһе last function сап be expressed 
in terms of a Fourier integral: 

eia . до 


деби = 5— [e TES, (3.11) 


the path of integration running from +00 to —c in the upper half-plane of the 
complex variable o. 


2 [n these formulae Vw? — wp? = |V w? — o,?| if © > од and Уо? — o, = Уи? — 0] 


if O < On 
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са 
t2 


Formulae (3.9) are linear in M; thus we can consider the fields Е? and Н? 
as a result of superposition of the fields Е"? and Н“ corresponding to steady-state 
waves caused by the radiation of dipoles with moments oscillating with various fre- 
quencies w and with various amplitudes; these amplitudes depend on œ and are 
determined by the fraction in the integrand of (3.11). After performing the superposi- 
tion we obtain 


Eo = Oe ee 3 hae к, um) Ф,» 


216 1 2 (3.122) 
AP = — — Qn UP Dm 
EC = = Мы Шу) АА ue) grad, Qm 
ات‎ U” n x grad, Ф 
E m tt 5 2 * m (3.12b) 
^ (m 10 я. 
HP СУ 225 U? grad, Фи 


1 
т=з TEES (P, ше — Аи? TOS uem n X grad, Pm 


where 


(m) еге 1 1 3 
== ---- em E ir) ойу ee ИИ # e 2 5 
и Кот Er рази в | ot A 0? — 0, | dw, (3.13a) | 
(т) eia 


Y | 
w 1 
Ш +--с | —exp—ilwt—7 yown] до, (3.13b 
Әлі "E ou WO р | . 2) 0), ) 
y 72) Om C + 


d lc sas dw, (3.13c) 


оре 
2ni | Vo*—o,3 ооу. 
L 
(m) 2A eia ih ч 
uz = sd = exp —1 C - 5. ү DEN do, (3.13d) 
i£ 5 


€ — 00 
(m) 0 eia ih А 
и = И 2zi ү 0?— Om —exp—i[ot—2 y o*—o,?| do (3.13e) 
ORAL w— 09 e < 
, L о Ы 
(т) А eia 1 ————— 
даа . Ў 2 
На i fe a exp — i (e t— - V а) do. (3.13f) 
5 | 


“Неге we have inserted z, = 0, which we shall assume throughout. We can always 
pass to the се of 207 0 substituting in all formulae = — z for 20. The function 1/0”) 
and its derivatives is obtained from (3.13 a) — (3.13 f) by replacing c, by 0,. 


Thus, the question is now reduced to th bl f i 1 
ШЕ INT e problem of effectively evaluating the 
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4. The function u™ and its derivatives 


We shall evaluate the five functions (3.13) following a method developed by 
Rubinowiez (1950, 1954). The path of integration L is deformed so as to make the 


exponent —i (o = DEDE) on the new path L’— purely immaginary. АП the 
с 


Imo L 


Fig. 3 


integrands in (3.13) have branch points +w,,. We place therefore a branch cut on the 
real axis of the w-plane between —w,, and +w,, (Fig. 3). The path of integration 
then proves to be an ellipse given by the equation 


Om 


= ————= (cos у + i В sin y), i 4.1 
о QR y T В еп (4.1) 
0x y < 2л, 
where 
Lg (4.2) 
в= 2 


et may be noticed that with a given constant value of z and with growing values of t this 
ellipse shrinks to the intercept (—«,, --о,) of the real axis. 
We denote now 


0° = o, t V1 В, (4.3) 
El eH 4.4 
= ЕВ , ( ) 
Е == уез, (4.5) 
(т) А 5m 
yo Sao. ES 2 Le | 
1 GA + о, : (4.6) 
yo” 


we then have on the‘ellipse (4.1) 


35 (= —2 o? — 2 = — in cos y, (4.7) 
C 
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„= (2+1), | (4.8) 

Vet ok = % (в 1]. 49) 

o ت‎ Е) (ғ- J ; (4.10) 
до =i ee € + dy. (4.11) 


We suppose that the pole 0); does not lie on the intercept (—w,,, + @,,) of the 
real axis, i.e. that there is not 


o = 0 and o < wm (4.12) 
Shrinking with the time the ellipse passes then the pole at «oy at the moment 
2 
{== сз 5 (4.13) 
v 


where 


бю ے‎ „Мощ — о*— On?) + (шї — 9? — On)? + 400520" 
V2 о 


It may also he verified (cf. Appendix) that at the moment t = 10”) there is 


y = ye. (4.15) 
Thus, after the moment 2 = tf”, a residuum of the respective integrand at the pole аж 
must be added to each of the integrals (3.13) when integrated not along L but along 
the ellipse (4.1). Let us consider the simplest of the functions (3.13), namely the 
function (3.13 d). We can write it down in the following form 


(4.14) 


2л 
ue = = Га ехр ( — 100 cos y)dy + rem (4.16) 
0 
where 
£2 К (т) И 1 
Am = ПН ры э л. у 
x + 5 д (т) ° (4.17 
70 Yo 
re” = Yt — 18) exp — би — = Va — es? — o). [R15 


We now denote 


та (m) | 
(m) D e yo д В 
fi 2 an TER NI exp (— и” cos y) dy ze тб A (4.19) 
0 
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2л 


eia 1 1 
fP = — exp (— i 2? cos y) dy. 
у” 


0 


Having in mind the known expression for the Bessel function of the n-th order 


2л 
Jn) = z- Г exp ( — im eos p + iny)dy, = (4.21) 
0 
we get 
u = oF Оре) +f + (4.22) 


To simplify the notation we shall omit below the upper index (m) in the symbols depend- 
ing on w,, such as f, fg, yo, (т, 10”, leaving this index understood. 

In order to evaluate the function f, we expand the fraction у (Е — yọ) in terms 
of thepowers of y,/é (if y > |Yol) or in terms of the powers of £/y, (if y < |yol). Denot- 


ing 
yore 13015 (4.23) 
where, according to (4.6) and (2.18), 


V — (02—02 — о„?) + Ve — o? — о?) + 400%, 


SP те (4.24) 
tg 0 Va 
= 5 «0x o) A 
tg 0— Е 
ГЕ сит ; (4.25) 
tg 0 + 0% 
we obtain 


— gia С — 1)" eind EN. 4.26 
л= У c [m 7.0), em 


n-1 
ja = — 1)" етіп E). % 
Ao ce Cre (Z) до 


+ exp —i (5, = = Vos — wm? — a) 1 (4.27) 


The first of these series converges better for y > |yol, the second one — better for 
y < lol; they both, however, converge for у 2 |yol. In a similar manner we obtain 
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an expression for the function f3: 
со 1 n 
fa = ей > (-і егін Li Jan). (4.28) 
SE y [yol 
This series is always convergent, as always у > 1, and |у > 1 (see the Appendix). 
We have thus explicitly evaluated the function u”, expressing it in terms of the 
convergent series f, and f. Before we pass to the remaining four functions (3.13) let 
us discuss two cases: 


ie 0 == O OOO 


In that case 


т) — E SET 4 4.29 
т = С ОКЕ b (4. ) 
0 


v denoting the group velocity of the mode in question. y, is a real positive number, 
and therefore 


0 = 0, [yel = Yo: (4.30) 
The residuum to be added at the moment NO = = becomes 
gr 
exp — i | Wot — = 002 — Om? — 2, (4.31) 
дч 0 == 0, Wo << Om 
We now have 
рез) (4.32) 
which is obvious, as the pole is never passed. Next, we have 
Wm 
tg 0 = — == l; [= (4.33) 
% 


In this case we can expand the fraction in (4.19) only in terms of y,/£, as now it is 
always у > |у. Therefore only the first expression for fj, namely (4.26) remains 
valid. 

Basing on the previous results we easily obtain the remaining functions (3.13). 
We have namely, in accordance with (4.8) and (4.9)? 


А 2л 
a eia 1 
uum rr с oe : A ein cos v dy 


Om 1 
2, ( -i) 


3 With regard to the sign of Иш? — те cf. footnote (?). 
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ІС 
a: C — 9). (4.34a) 
D -— 77 
2n 
Or 1 
eia D (: xis i 
uem ES ec 2 А ein cos y аф 
2л о), 1 
са 


та А 1 ( 
сете dum 111. - Wo m eeu егіп cos v dy 
an Wy? — Wn" 70 ув ее 


Yo 


Op 
= س‎ с | J, (0) + = dale (4.341) 
с | e Jo (1) terra fa 
2л 
Ou а On PE И Кез о, Yo 
e 2л Фе F) "E aeo £— yo 
1 1 —in cos St Та iamnm 1 
старо rdy =i |e oF (+) ле 
C 
IE e AA aM roy ea (4.340) 
2л 
и = = | А ene v dy = ei Л (M) + А + fo (4.344) 
0 
Qu 2л 
(т) __ дена От р ей От 


1 А -incosy dy = — —— © кара 
ue CG 20 2 Dre s PTS 2n А 2 Е 
0 0 
== Vo? — о? "ERA CONI SU UN E dn етіп cos y аф -— 
712 
(, " 1) ло) + Vo? — On? (А -9| (4.346) 
À 


g 


ги _ Om 
fs л) + esr 
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2л 


та 1 А 
ip =i [= (ret aenmas = 
0 
2л 
eia w 1 ^ А 
acm EE ES pant e A -in cos 
== IE (: Е) + ay] vdwy 
0 


-—i E си. | y— D (n) + её do Jo (0) + do (A+ | . (4.349 


In order to obtain the field of the m-th mode, we must шеп ше. теа] рагїз оЁ 
(4.34) into (3.12a) and AG 1р). Putting о = 0, « = 0 we obtain po and H т). Putting 
e #0, «= M the field E". H qoi is obtained. Finally, inserting о 72 0, х = т/2 we obtain 


pm and H (m), In turn, the field E™, НХ? can be found according to (2.15). Sum- 


res ? res 
ming with respect to т yields the final result. 
We can now find the field radiated inside the wave guide by a thin antenna of 
finite length supposing the current in the antenna to be 


Кхо yo, zo) * (t) B еге sin ct + (1 — e?) cos ont | 
o n 


(cf. 2.14), where xy, ys, 2 are coordinates of a point on the wire. 

ү, : I (хо yo» 20) 48 : МЫ Re : 

e replace M in (3.9) by mu A Thefields E, H, E, H, E, H as obtained 
a 

from (3.8) and (3.10) are then integrated with respect to ds,. If the coordinate z, 
varies along the antenna, z must be replaced by z — 2, in all formulae before 
performing the integration along the wire. 

At distances far from the wave front the need of asymptotic formulae arises. 
We shall derive such formulae in the second part of this paper. 


Appendix 
We shall prove that with y, defined by 


o © 
Yo = sa sls да --1, | (А.1) 
we have always lol > 
Let us denote 
_ De _ %-і% 
= or d (A.2) 


(A.1) then becomes. 


je 5572 a Vz? — 1. | (A.3) 
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This is a transformation transforming the exterior of the intercept (—1, +1) of the 
real axis of the z-plane into the exterior or interior of a circle of radius 1 in the yọ- 
plane. According to the convention concerning the sign of the quadratic root which 
was adopted in footnote (2), we choose the first of these two possibilities (Fig. 4). The 
intercept (—1, +1) itself is transformed into the circle of radius 1. Denoting 


Fig. 4 
yo = х + iy = |у) е, (A.4) 
г.и |- по, (А.5) 
it is easily seen that 
1 1 
=> пе 50; А. 
и 9 (њ чс я cos 0; ( 6a) 
1 1 
== — — | sin 0. A.6b 
E i 
Thus an ellipse with semi-axes 
1 1 H 1 ) 
= — mss b= — — — A.7 


is transformed into a circle |yọ| = const. The points A’, В, С’, D’ in the yo-plane 
are representations of А, B, C and D in the z-plane. 

Now, z can be either real (о = 0), i. e. lie in points like B (w> о) or P, (о<0),)), 
or complex, i. e. lie in the lower right quarter of the z-plane, e. g. in a point like P, 
(020). B . . 

We see that in the case of 0 = 0 and w,<w,, there is |yọ| = 1 (point Pj"); in all 


other cases there is |yy|> 1, which proves the thesis. 
Next, we derive equations (4.24) and (4.25). We multiply (A. ба) by sin 0, and 


(A. 6b) by cos 0. Adding and subtracting the results we obtain 
u sin Ө = |yg| cos 0 зад, = (A. 8a) 


1 
v cos Ө = — cos 0 sin б. (А. 8b) 
lol 
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eue ge PE аа ЖАУ ЗА RE НИ ес. 
Multiplying (A. 8a) by (A. 8b) we obtain an equation from which (4.24) follows imme- 
diately. Eq. (4.25) is obtained by dividing (A. 6b) by (А. ба) and solving for |yol- 


2 LA 
We conclude by E that at the moment #7 = ход NC there is 


y = Inl. | (4.0) 
(сї. (4.15)). 


At the moment (” we have 


дэ лы | 
кые S (4.10) 


where 
2 

Во 1 Оп) x 
From equations (A. 6) we have E 
u? v2 - с | 
пет Це еы т 1; | 
: Bu (All) А 

+ is : 
4 (bo m a т (bi Po 2) : 3 
At Ше moment t = 409 the ellipse (4.1) just passes the pole oy — ig and thus eq. | | 


(4.1) yields 


Wy — 10 = 


—— (cos y + ша y» - HA) 
А | g Е ion р 
and it follows that өтт үш a= al i M 


cur i iie mm ле 
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КРАТКОЕ СОДЕРЖАНИЕ 
Гаевский, Радиация диполя в волноводе в неустановившемся состоянии I 


Поданы точные формулы, определяющие электромагнитное поле, происходя- 
зцее от диполя, помещенного в волноводе любого сечения, если ток в диполе 


зависит от времени, как 1 (1) (1-е г) ты 


И REFERENCES 


de Broglie, L., Problèmes de propagations guidées des ondes électromagnétiques, Paris, 1941. 


Cf. also Russian translation, Moscow, 1948. 
Cotte, M., C. R. Acad. Sci. [Paris], 221, 538 (1945). 


Chien, W. Z., Infeld, L., Pounder, J. R., Stevenson, A, F., and Synge, J. L., Canad. J. 


Research A, 27, 69 (1949). 
| Cohn С. I., Proc. Nat. Electronic Conference, 8, 284 (1953). 
Gajewski R., Bull. Acad. Polon. Sci. IV, 3, 29 (1955). 
` Kahan T. and Colombo S., C. К. Acad. Sci. [Paris], 226, 2060 ы 
Poincelot P.; Ann. Télécomm. 9, 315 (1954). 
Rubinowicz A., Acta Phys. Polon. 10, 79 (1950); 13, 115 (195%; 14 (1955); 


Zeitschr. für Angewandte Math. и Phys. in Ше press. 


= 
me € 
er pat rents 4 ы | Е 
Varios b iT 8 р, 22 te "i " LIN, 
2 Jud. xk аа "cx و‎ же = f 
ж he Cor dil er LE | масс 


ЖЕ дано " DIM 27,03! FT. j 
E BE st зале PES 


eo وف‎ " 


Vol. XV (1956) ACTA PHYSICA POLONICA Газе 


ON THE LONGITUDINAL ATTRACTION AND REPULSION 
OF SPHERES IN VIBRATING AIR 
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The force acting on spheres in a Kundt tube with their line of centres parallel to the 
tube axis has been investigated. As in the previous paper on the mutual influence of 
spheres in the transverse position it is shown that for the longitudinal position this force 
depends on the frequency and intensity of the air vibrations as well as on the size of the 
spheres. The force at very small distances is shown in some cases to be one of repulsion, 
so that Andrade's results cannot be generally true. 


The explanation of the dust tube phenomena in the Kundt tube given by Kónig 
(1891) was shown by Andrade (1931, 1932) to be inacceptable on several grounds. 
Andrade showed that the viscosity of air is essential for the formation of ridges and 
studied the forces acting between cylinders or spheres in strong acoustic fields. In 
the course of his investigations he arrived at results that contradicted in part the former 
experiments of Cook (1902), viz., he found that two cylinders or spheres situated trans- 
versely to the vibration vector attracted one another at distances smaller than their 
radius, whereas according to Cook they repelled one another at such distances. It has 


been shown by Staszewski (1954) that the forces in question are largely dependent on 
` the size of the spheres аз well as on the frequency and intensity of air vibrations, so 


that the contradiction of Cook's and Andrade's results may be simply Mies by 
different conditions of their experiments. 

As concerns the forces between spheres situated parallelly to the vibration vector 
Cook and Andrade agree that there is attraction when the spheres are very near and 
repulsion at larger distances. According to Andrade at still larger distance, at which 
the vortex systems around the spheres just touch, the repulsion turns again into weak 
attraction. At small distances the longitudinal position becomes highly unstable the 
spheres tending to set themselves transversely to the vibration vector. 

The attraction of spheres in longitudinal position takes place at distances smaller 
than the radius (Cook), that of cylinders at distances comparable or less 
than their diameter (Andrade). Neither of these authors studied this situation more 
closely. Our experiments on the mutual influence of spheres in the transversal position 
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having shown that the force at very small distances may be either one of attraction 

or of repulsion it was thought worth while investigating also the behaviour of spheres 
in the longitudinal position. 

A study of the dependence of the force acting between spheres with the line of 

3 == centres parallel to the vibration vector 

| was made using the apparatus descri- 

bed in the previous paper (Staszewski 


1954). 


2 
Fig. 1 Fig. 2 


Fig. 1. Displacement of the suspended sphere from the dynamical zero ZZ. 


Fig. 2. The suspended sphere with the stabilizing device. 


А sphere 5, was suspended in the Kundt tube about 4 cm wide near Из open 
end on a fine glass fibre (0.03—0.04 mm). Another sphere 5, was sealed to a thicker 
glass thread attached to the arm of a micromanipulator and could be brought to any 
position with respect to S,. The positions of the spheres were determined by means 
of a measuring microscope of a Zeiss comparator for spectrographic work. The displa- 
cement of the sphere 5, was measured from the dynamical zero, i. e. from its position 
of rest ZZ (Fig. 1) when the vibrations obtained in the tube, the sphere S, being re- 
moved. It served as a fair measure of the force between the spheres. 

On the ground of symmetry the force acting between the spheres in the longitudi- 
nal position must be directed along the line of centres. If this line deviates somewhat 
from the vibration vector another component, perpendicular to the line of centres, 
appears at small distances and the sphere 5, moves around the fixed sphere S, until 
their line of centres becomes perpendicular to the axis of the tube. In order to avoid 
this instability the movable sphere was constrained to move in the direction parallel 
to the axis of the tube by means of a transverse fine glass fibre F, (Fig. 2) attached 
to the brass ring supporting the glass fibre F on which 5, was suspended. The stabiliza- 
tion in the longitudinal direction was obtained by immersing the glass fibre F, in 
a mixture of honey and glycerine. The stabilization for the smallest spheres used was 
rather bad and many (up to 20) readings had to be made at each position ofthe fixed 


sphere. The average displacements of the suspended sphere gave, however, rather 
smooth curves. | 
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Let us consider first the influence of the size of the spheres on the force between 
them in the same soundfield at the frequency v = 360 c.p.s. and at 60V potential 
difference on the voice coil of the loudspeaker. The spheres had diameters from 
0.6 mm to 2.6 mm. The results of the measurements are given іп the graphs of Fig. 3. 
The force of repulsion is considered here as in other instances as positive, that of 
attraction as negative. For larger spheres (uper curve) the repulsion increases con- 
siderably with decreasing distance of the spheres up to their contact. For smaller ones 
the repulsion increases with decreasing distance down to about 0.5 mm and then 
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| attraction 7] 2 3 тт distance 


Fig. 3. Displacement of the suspended sphere S, as a function of distance from the fixed sphere S, for 
spheres of different diameters at constant frequency v — 360 c.p.s. and constant voltage U — 60V 
on the voice coil of the loudspeaker. 


1 2 3 4 mm 


Fig. 4a Influence of the intensity of air vibrations on the force acting between spheres of diameters 
d, = 3.1 mm, d, = 3.2 mm. 


varies but little up to €ontact. The repulsion of the smallest spheres increases with 
decreasing distance, passes through a flat maximum, falling to zero and changing into 
attraction when the spheres nearly touch. This behaviour is opposite to that in the 


46 В. Adamczyk and W. Staszewski 


Кы COCA Feuer. e = Е = 


transverse position, where under similar conditions the smallest spheres repel at con- 
tact and attract one another at larger distances (Staszewski, 1. c. p. 222, Figs. 20 and 21). 

It is to be noted that the forces represented by the curves of Fig. 3 are not directly 
comparable inter se, the coefficient of proportionality between the On and the 
displacement being different in the three cases owing to different weights of the 
spheres, since for heavier spheres larger forces are necessary than for the light ones 
to effect the same displacement. Corrected for the weight of the sphere the uppermost 
curve would rise more steeply with decreasing distance of the spheres. 

The influence of the intensity of the acoustic field is shown in Fig. 4a, 4b, 4c. 


mm 
0,4 


mm +. (/- 105 V 
+. U= 105 V 


о—о—о—о= 60 V 0,2 


о-оо 60 V 


0,2 repulsion 


1 attraction 2 


Fig. 4b Fig. 4c 
Fig. 4b. Influence of the intensity of air vibrations on the force acting between spheres of diameters 
d, = 1.5 mm = d,. 
Fig. 4c. Influence of the intensity of air vibrations on the force acting between spheres of diameter 
d, = 0.6 mm = dg. 


For larger spheres (Fig. 4a) of diameters d, = 3.1 mm and d, = 3.2 mm keeping 
the frequency » = 360 c.p.s. constant and raising the potential difference on the 
voice coil from U = 60V to U = 105V entails fairly uniform increasing of the repulsion 
at different distances. 
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<2--е--е--%- /--360 c.p.s. 
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1 2 3 mm ` 


Fig. 5a. Influence of the frequency of air vibrations on the force acting between spheres having the 
diameters d, = 2.5 mm; d, = 2.6 mm; frequencies v = 360 c.p.s. and ? = 530 c.p.s. 


In the case of smaller spheres (d, — d, — 1.5 mm) the repulsion attains the 
maximum at about 0.5 mm distance and decreases uniformly when the spheres are 
brought still nearer (Fig. 4b). 

For the smallest spheres (d, = d, = 0.6 mm) the maximum of repulsion occurs 
at the larger distance of 1.6 mm (Fig. 4c). With decreasing distance the repulsion 
decreases to zero at about 0.5 mm and turns into attraction which ‘increases steadily 
up to contact of the spheres. 


The attraction of very small glass or wax spheres (0. 2—0. 3 mm) in the Kundt 


ы теле © 
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tube was observed by Cook at distances smaller than half their diameter, but it is 
not clear from his paper whether larger spheres showed a similar effect. Andrade found 
that if two cylinders of 1.5 mm diameter were very near together the suspended 
cylinder tilted towards the fixed cylinder before it was displaced sideways, showing 
that there is an attraction at small distances. However, in view of our measurements 
Andrade's statement that in the longitudinal position *when the obstacles are very 
near, there is always an attraction* cannot be generally true. 


pulsion 


1 


20.2 attraction 


Fig. 5b. Influence of the frequency of air vibrations on the force acting between spheres having the 
diameters d, = 0.6 mm = 4;; frequencies v = 360 c.p.s. and v = 165 c.p.s. 


The influence of the frequency is apparent from Fig. 5a and 5b. Fig. 5a gives 
the force acting between the spheres having diameters d, = 2.5 mm and d, = 2.6 mm 
as a function of distance for the frequencies у = 360 c.p.s. and v = 165 c.p.s. at 
constant voltage U — 60V. Fig. 5b gives this force for small spheres (d, — d4 — 0.6 mm) 
for the frequencies у = 360 c.p.s. and v = 165 c.p.s. 

For large spheres the repulsion at contact and, for the small ones the attraction 
decreases with increasing frequency. This effect is probably largely due to the larger 
amplitude of the vibration of lower frequency when the voltage on the voice coil is 
kept constant. 

It would be worth while investigating the influence of the amplitude of air vibra- 
tions measured directly in the Kundt tube by the smoke particle method of Andrade 
and Parker (1937) which is, however, not an easy task at least for higher frequencies. 


КРАТКОЕ СОДЕРЖАНИЕ 


Адамчик и Сташевский. О продольном притяжении и отталкивании шариков в вибри- 


‚ рующем воздухе. 


Были исследованы силы действующие между шариками помещенными в трубке Кундта 
при параллельном положении линии центров по отношению к оси трубки. Обнаружено, что 
сила между шариками может быть равно притягательная, как и отталкивающая в за- 
висимости от частоты и напряженности звукового поля, а также от диаметров и расстояния 
шариков. Обнаружено, что в некоторых случаях отталкивание бывает и при очень малом 
расстоянии между шариками, вопреки утверждению Андраде. 
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GENER ALIZATION OF THE METHOD OF SUPPLEMENTARY 
VARIABLES TO SYSTEMS COMPOSED OF TWO 
KINDS OF PARTICLES 


Ву Z. GALASIEWICZ 
Institute of the Theoretical Physics of the Polish Academy of Sciences, Wroclaw 


( Received. July 30, 1955) 


It is shown that the elementary excitations in a system composed of electrons and 
ions (interacting electrostatically) are of two types: a boson type and a fermion one. The 
boson type excitation may be described as a system of oscillators vibrating with the 


frequency үш? + we, where оу is the Langmuir frequency of the electron vibrations 
(vibrations of the electron plasma) and с), is the Langmuir frequency of the ion vibrations 
(vibrations of the ion plasma). If one neglects the (Coulomb) interaction between the 
ions and electrons, the elementary boson-type excitations are given by two systems of 


oscillators vibrating with the frequencies w; and w. As V ws + ws < 0 + 05, the energy 
of the zero-point oscillations in the case of non-negligible interaction between the ions 
and electrons is smaller than when this interaction is absent. The energy difference of 
both these zero-point oscillations may be considered as a measure of the binding energy 
of the system coming from the *long-range^ Coulomb interaction (which is responsible 
for the boson-type excitations). This energy has been evaluated for Li, Na, and K and 
compared with the total binding energy calculated by other methods or found experiment- 
ally. The orders of magnitude are the same and the results of the author lie nearer to the 
experimental values as those following from the Hartree-Fock method. 


1. Interaction 


In recent years many papers appeared concerning the multielectron treatment 
of the plasma and the solid state. These investigations developed especially in the 
USSR (Vlasov, Landau, Vonsovskii, Bogolubov, Zubariev, see the Review by Bonch- 
Bruievich 1955) and recently also in other countries (Готопара, Bohm, Gross, Pines 
et al., see the references in the papers of Ecker 1955 a, b). The most promising method 
seems to be that proposed in two different forms by Zubariev (1953 a,b) and Bohm & 
Pines (1953). . 

In the present paper Zubariev's method is generalized so as to be applicable to 
systems consisting of two kinds of particles: electrons (fermions) and ions (fermions 
or bosons). In another paper of the author (to be published in this journal) a similar 
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generalization has been performed for the BP method. The results obtained by both 
methods turned out to be coincident. It has been also proved that both methods 
though not identical, are strictly equivalent. № may by noted, however, that they 
offer somewhat different computational possibilities and also different possibilities 
of interpretation. 

The main idea of Zubariev’s method consists in introducing into the electron 
wave function of supplementary variables gy, which are Fourier components of the 
density operator. The wave function depends on the 0, components with indices k, 
for which k =| k| < ky. This cut-off allows the separation of the collective part ofthe 
plasma movements. The quantity ky must be finite but otherwise it remains for the 
time being undefined. In consequence of the equivalence proved by the author the 
value of k, is here taken equal to that of k, used by Bohm and Pines (1953). The 
considerations leading to the evaluation of k, apply then also to ky and can be easily 
generalized to systems consisting of electrons and ions. After introducing the supple- 
mentary variables the Hamiltonian can be split into three parts refering respectively 
to the energy of the fermion-type particles, the energy of collective vibration of the 
boson type and the energy of interaction of the particles. In the case of Coulomb 
interaction the second part of the Hamiltonian represents oscillators vibrating with 
the so-called Langmuir frequency (Langmuir and Tonks 1929). In Zubariev's paper 
(1953a) the Hamiltonian was only approximately hermitian. Migdal and Calitski 
(see Zubariev 1953b) modified the method of supplementary variables in such 
a way as to obtain a strictly hermitian Hamiltonian. This method has been applied 
in the present paper. 

All the calculations are performed in the approximation in which о), w, and œg 
do not depend on X. This corresponds to omitting in the Hamiltonian the term describ- 
ing the interaction of fermions with the collective motion (of bosons). The taking 
into account of this interaction would require a special investigation. 


2. The introduction of supplementary variables into the Hamiltonian 


We consider a system of interacting particles consisting of N, electrons and М, 
ions. The Hamiltonian of this system has the form 


- 2, SS ON ETA Мм 


АЖА . ll, (1) 


1 
+ D W (Fu — в) + Y Us — Ну). 
lh 7,1 
Неге ғ» p; and m denote respectively the coordinates, momentum and mass 
of the electron, R, x, and M the corresponding quantities for the ions, N, is the 
density of the electrons, N, that of the ions. The three double sums represent successi- 
vely the energies of the electron-electron, ion-ion, and electron-ion interactions. 
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In the method considered the wave function depends оп г, о, R, X (coordinates 
and spins of electrons and ions) and on the supplementary (redundant) variables 
0, (Е < ko) describing the collective motion of the particles. The wave function has 
to satisfy the subsidiary condition (see Zubariev 1953b) 


Op lin cer e na Seer ML OSEE (2) 


We shall show now that thanks to condition (2) the Hamiltonian (1) can be given 
a form in which appear not only the variables Qr but also their canonical conjugates. 


Thus we shall single out from the Hamiltonian the part referring to collective mo- 
tion. | 


To this end we expand the potential energies into Fourier series, as follows 


И, — n) = Хоф 
k З 


W(R, — R,) = 2 w (k) е (Fı, Pi) : (3) 
Ullr; — Rj) = Mu (k) e" 07—80 
k 


For Coulomb interaction, with which we shall exclusively have to deal in the 
sequel, the coefficients of these expansions have the form 


4ле? AnZ?e? . 4л2е?2 


v (Е) = Eye w (k) cu k2 , u(k) T Ta (4) 


where —e is the electron charge and Ze the ion charge. We introduce the density-of- 
charge operator 


М, № | 
Ona Mie е Аы NE of с (5) 
т 1/1 k 


By E the corresponding components of the expansions, we get 
Ni NUE 
Qum D dap Ze Ме ори ооа (6) 
ІЛ JA E 
When the system as a whole is neutral N, — ZN, and the mean density of 
charge is : 
о = -еМ + eN,Z = 0 : (7) 
We transform Schrüdinger's equation ih 99/9t = Н, and condition (2) by 
means of the unitary transformation 


i 9 
U = Tr, Ex eae , where Р, = — ih 20, (8) 


and о, is given by (6). The subsidiary condition may be written as follows 


О.ф = 0, | ке) 
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where 
0, = 00,0%, Фф = Up 
If we take into account the relation 


U+ (Pa) дк — Qu Ut (РӘ = — ih 55 ОР), 


we see that 
, , - а P А 42 p 0 , , 
Qu еже d ге cd Ut e = (Qr — e) 9. 
Thus under the unitary transformation (8) the subsidiary condition (2) becomes 
(Qu — 00% = (О, — о? — оо = 0 for k< ko (10) 


We are going now to subject our Hamiltonian to the transformation (8), making 


use of the following relation 


p; Ке) = fr) pj — 2ih (Vs; fr). р) — ^? vs fir). (11) 
then 
М, 1 
сай! Y Up? Ut = = р? те Е hk 
2m 2 Pj А йй шй Е 
1 ill ГА kk, 
са (12) 
паа Bi Pr, (ЁЁ) е-* К») a7 
2 oss 
kak: < ky 


The analogous terms relating to the ions are obtained by the substitution 


(bj. En -e) ED. (В), л, М, Ze) (13) 
Hence 
1 N: М, 7 5 à T 
Sas ре а 7t; LZe > 
DS Un; U = 2M vd № е mn ( =) 4- 
a JE Lk, ke hs 
5 ом | (192) 
d L kı, ki< ko Р.Р. (А, ko) e^ i (ch) R; 
ka, k< ko 


We split the Coulomb interaction in “long range“ (k< ko) and “short range“ 
interactions (k> ko) 


1 
5 У? V (|r; — rl) а Da W( в, — В.) + 23 U (|r; — в) 


АА L+h ET 


UEXESeeayytzee. 


Jit k, R< Ro Tij. k, kh. 
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1 2 2 
= а емінен) 2 ee eik (tj, Fı) 

lh К, k< ko LEE k, k>ko 
AnZe? | 2 

E. № У ща ек (rj—Rj) = 3 3 tu е, (;- Ry) 

j,k k, R< Ro 251 k,k>k, k 
v(0 w(0 
+ 0 e+ UD) mom (14) 


Аз w(0) = 2%(0), u(0) = - 70), and №, = ZN, (system neutral as a whole), 


the last three terms cancel each other: 


0 0 
— м m Z Ni — «0Z2NN,—0 (15) 


Similarly we can split the charge density in d eu and “short-range“ densities: 
e(r) = 2 Фе EDD 2 
Е, kk. 


It follows from (6) and (10) m 


On P = 0,% - (e eis +2 еви) gy’ for ЕФЕ 


jA 
In (14) the terms relating to the long-range“ Coulomb interaction (в< ko) can 
be expressed as functions of the Fourier components of the density operator 0, 
(k < ko). In consequence of condition (10) the supplementary variables Q, describing 
together with the variables Р, the collective motion (collective vibrations) аге intro- 
duced into the Hamiltonian acting on the wave function. 
Taking account of (6), we calculate the expression 


1 4л 1 4ле? ^ 
а — —— کے‎ (5 715) 
2 е8 0 2 № pei 
k, k< Ко Ji; Je К, k& Ко 
1 47? е? " AmZe? . 
— ксы ik (RyRy) = ik (r;—Rj) 
рУ гре еі es TO pres ve (16) 
1, ВЕ, kx Ro 1,7 К, Е< № 


which differs from the “long-range“ part of formula (14) in the terms with д =} 
and /, = 4. These terms have the form 
N, 
M 2 
2ле? У 72 + 2nZ?e? 3 gi (17) 
k, k< Ro k, R< ko 
On subtracting them from expression (16), we get the “long-range“ Coulomb 
interaction in terms of the Fourier components of the density operator: 


1 М N, 
2л т Ок -k — 2ле? M ii — 2nZ*e? 3 72 3 (18) 
k, k< № k, k< ko k, k& ko 


у> 
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Thanks to (12), (12^), (14) and (18), we see that the Hamiltonian (1), after applying 
to it the transformation (8) may be stated in the form 


71 1/1 jı Fis k, k> Ко 
1 |. 4nZ?e? ik (Rj —R}) 4nZe? | 
l*Fl; k,k>ko Jl k,k>ko 


ем , Z?e?N, 1 
[Sm ый, D, Рада У, erent 


К, k< ko k, k< Ro 
— 2ле? у та гл ?е? 1 72 d 
k,kR<k ” К, k& Ro 
ie б hk iZe ; hk 
RE етік" р ES P: eikR P Rips 7L + 
m ; | Jik 2 J ae M ) 1Ғк 2 1 
j,k, k< ko l, k, k< ko 
s PIRE icm 1226 3 Р, Ру, (Е By) ei a tho R 
— enter the) тр — + е за та) 287 
Әт № kı к, ( 1 2) J 2M ktk 1°2 
ЛВ, xh lu kis Fı < ko 
a» Ка & Ко ka, ke< Ко 
ki+ki +0 кк, +0 


The summation in the last terms of formulae (12) and (12’) has been decomposed 
into two parts: for №, = —k, and Ау Æ — Бо. 
3. Further transformations and. discussions of the Hamiltonians : 


The terms in the Hamiltonian not depending on the coordinates of the particles 
have the form 


e?N,- ZEN, 1 
Бе Е S x] à, БАРА? У genes (20) 


k, ESko k, k< ko 


and after taking condition (10) into account they become 


H N Z3e* N, 2л 
е 1. 2 9 
I 2 ( 2 2M | L P, D k2 Ок О (21) 


Оп changing to new variables 


А | | 
Һ- аттар» Qu= V2 А, qm (22) 
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where A, is defined by the condition 


e NIE a Lice NAT 
?m 2M | 2 = т^ (23) 
we get finally 
E Ame?N, | A4mZ?e*N, 
Нос = p \3 у А, (24) 
К, kS Ry 


As 4леё№ [т = оу is the Langmuir frequency of the electron vibrations and 
4nZ?e?N,! M = а» the Langmuir frequency of the ion vibrations, we have 


2 2,2 
of a RAR Vol tal (25) 


m M 
If N, = ZN, = М and we introduce the reduced mass u by 
|| T Z 
hat «аам Las 26 
"PIE (26) 
we get 


It follows from (26) that u = m/(1 + Zm/M) == m, as 2 < Мт for all physicall 
possible cases. 
We pass now to the boson number 
occupation representation: 
1 + i ^ 
== (hb b8), = —— (by — 6 . (28) 
дк үз (br + 82), Pr үз (b. — br) 
Then the operators b, satisfy the following conditions 
by, bj, — by Ôm = Oa 
b, bp, — by би, O (29) 
bg bt — bt bt = 0 
and . 


1 

HS = ha’ 5 (oi by, + 3 . (30) 
k, k< ko : 

This is the Hamiltonian of harmonic oscillators with frequency о). They describe 


the collective motion (collective vibrations). 


The term 
f М, 2 № 2 
- Р; M (31) 
set gute 
Ват m 
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where 0 1 4ле? 
Hi. = E N es: Gite a 
J jı ts > ko 
w(0 i Anze? . 
n S мы), Y екст) + (31) 


14-1, k, R> № 


u (0) N, Ny — У, pm TEE oM m 


jl k,k>k, 


(see (15)) describes the energy of a system of particles possessing the same masses and 
spins as the particles initially considered, but such a charge that the potential energy 
of their interaction is the “short-range“ Coulomb interaction. 

The interaction energy of the particles with the collective vibrations is described 
by the following terms 


H, = HO + Н® + HO) + HOO 


1e hk iZe hk 
niu pe ТН. (% Б 2 и ee (2 - =) 


ГА k< ko LE, k< ko 


22е? е 
E = № Pp, Py, (k №) ево"; — 3M 3 Py, Pp, (By k) ei th) Ry 


Шы ere (52) 
ki+k: +0 ki+k:+0 


4. The perturbation theory 


We shall deal now with the lowest energy state, which corresponds in the case 
of collective motion to the ground state of the oscillators. 


Using the perturbation formalism we have to compute the following matrix 
elements 


(АН), (АН), (АН?) (АН) (33) 
To this end, — similarly to what has been done with Hos, — we shall substitute 
new variables into (32) making use first of (22) and then of (28). 
(3) ehr E r: 
gi 5 mo’ X (ек, Ex,) (by, br, ср 55, bi, bn, by, тр br, b_k,)e puo Ч 


J, Ey RS Ro 
2» „< Ко 
ki+ka+0 (34) 
€, 13 a unit vector pointing in the direction of k. We shall denote the energy of the 


initial state by E, and that of the intermediate state by Err, then 


р? ресе ВЕ Fk) 
Е-е: Er, r, = 2h0 са iiu 
E 


(35) 
к — Eg ^ 280” 
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The corresponding matrix element will have the form 


Pme – у fas к= 
Ка, ki< ko AS 20) 
Ra; Ка< ko 


The matrix elements describing transitions to intermediate states are, by (34), 


2 
(3) ehr В $ 
(Hî ) оожак, = ; En, Ep, € ta) ty 
то 


(37) 
н e?hz м 
( T )к,к, :00 = 7 Ek, Ek, e-i(kı +k») ту 


, 


Substituting (37) into (36), we get 


ehi 2 2 1 , 3 у 
(ЛН) ооо = — es N, D fat — 18 (s) E EU deny 
кав 


Ву similar calculation we obtain 


, Е 
О ere [m В Дш (39) 
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where п’ is the number of terms in the sum У) and may be expressed in terms 


k, k< ko 
of ky by the formula 
mdr MA me 
SET (à) Done (40) 


In consequence of оу = а”, @ w’ and (as we shall evaluate presently) 
п’ == N (for М = М, = №; Z = 1) the terms (38) and (39) may be neglected in compa- 
rison with the zero-point energy n'hw’/2. 

We shall now turn to the remaining terms (32). As before, making first use of 
(22) and then of (28), we get 


2xh\ ^ hk 
apt y (2) СЕЛІ 


т + 
Js k, k< ko E (41) 
17 
72 ДЕ 4 
ЕШ ЕНЕ 


The energy of the initial state beeing Ер = р22т, and that of the intermediate 
state E, = hw’ + (p; — hk)?/2m, we have Ep — Ер = Во 
The matrix element looked for has Ше form 


+ 


Us Б, k< № 


f HO) (Н), 
(АН оо = — » MES EE ae , (52) 


. Jk, RK Ro 
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The matrix elements describing transitions to intermediate states are, by (41) 


» hk\ . 
(Ноль = = [a | Ek (» = м ei 


т 


(43) 

ev] oh A ES hh 

(Hf )ujo e | "2 | вие "е Ç x w), 
Substituting (43) into (42), we get 
|48-%21 
2ле?ћ СА ООЙ 
(ЛН) оо = — аня 

т?ау , 2 " ho (44) 


дата (m) ту RER" 
TORNA " 2m 40 т у 


By similar calculation we obtain 


(2) natos É л? 9 as) 

Е a EJ D (%- 0 M | 9) 

Аз 05 < w’, the term (45) is small in comparison with the kinetic energy of the ions. 

It follows from our calculations that the most essential contribution to the lowest 
energy state comes from H$”. Bohm and Pines (1953) showed that this term can be 
eliminated by means of a suitable unitary transformation; the new Hamiltonian refers 
to particles with a greater effective mass and the frequency w’ becomes a function 
of k. If we evaluate the order of magnitude of п’ (or k, = ko, see the Introduction) 
basing on the requirement that the energy of the lowest state be as small as possible 


and on the assumption that М, = № = М, Z = 1 (similarly as it has been done 
by Bohm and Pines (1953) $ 2) we get п’ == М. 


5. The case of noninteracting electrons and ions 


We shall consider now a system of electrons and ions which do not interact with 
one another. The Hamiltonian of such a system has the form 


М, р? ек (rj, A) Y , 1 
> 2 + 2ле? 3 TN 2л Ме? po + 
jn Ран! E |. (46) 
№ 2 А 
It] г ек (Е-Е) С 
epi n PD у, три 
11 1+ 1,6 k 


The apostrophe is used to indicate that the zero term in the expansion of the 
Coulomb energy has been dropped. This corresponds to the fact that we are con- 


sidering the interacting electrons on a background of the neutralizing charge of the 
ions and vice versa. 


EA 
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— M ÀÀÀ——. 


For the system described by the Hamiltonian (16), we are looking for a wave 
function of the form 


Ф = Pa (91... гм, ом,›1/...00-..) Pa (Ry Zi. Ry, Ext]. QD). (47) 


The supplementary variables Q® describe the collective motion in the “electron 


system" and the variables Qf? in the ion system. 


If we assume that О, = Qj? + 02 (see (6) and (10)), the subsidiary condi- 
tion (2) takes the form 


Qn p = (00 + 09) Pa Фа (48) 
and splits into the two following conditions 
Qv —-0; QP o = 0 (49) 


We introduce now the density-of-charge operators for the “electron system“ and 
the “ion system“: 


Ni 
00) (г) = — e 3} д (r — г) + 2№,е; o? (r) = Ze а ó(r — R) — Ме. (50) 


jll 


The corresponding Fourier coefficients are 


NS АГ НЕМИ: É 

Прат teh вае САА Ты A, 2 ЕДЕН (51) 
НА 1 | 
00 =0 ; 00 =0. (51) 


We transform next the Hamiltonian (46) and the supplementary conditions (49) 
by means of the unitary transformation 


2 i PY 00 4 р (002), 
е, да. KU OR EAE RA DE (52) 


where 


PË ih 50; (QU) as за 1.2) 


The form of the transformation Ü shows that both the parts of the Hamiltonian 
(46) transform independently. The subsidiary conditions (49) become 


(QP — oP) =0 ; (00 – 00) р = 0; EH (53) 

where | 
y = Пр. G=12) (53^) 
We shall express the “long-range part“ of the Coulomb interaction between the 


electrons in terms of the 00) components, and the corresponding part for the ions 
in terms of Ше 000 components, similarly as in equation (18). Then, thanks to con- 
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ditions (53), we shall be able to single out from the Hamiltonian (46) transformed 
by (52) the part 


HM = un У G PË po (D 25 op 04) 


2m 
k, kS ko (54) 
2,2 
1 2 Ы n | k2 PË p + (2) a 
E k, k< ho 


We introduce the new variables 


h i i а . 
Ро = ра ; e = V2 PP qi, (2 Les 1,2) (55) 


V2 29 
where ; are defined by 
2 262 
LA 1802 = a ал (aus ; Ea Дара = 3 (ДЇ). (56) 
Then Н’,, becomes 


ho , (1) pa ау (D ho , (2) (2) (2) (2 
H ose = 2 E (pk p-k-- qk 9-ю) + —у (Pk P-k + qk 4-4) (57) 


where w, and а»; аге the Langmuir ее of the Cu and ion vibrations, 
appearing already in formula (25) for о’. 

By means of similar relations to (28) we pass to the boson occupation number 
representations of two types, in which Н, will be expressed by 


: 1 1 | 
H „= (D+ D) Er (2)+ 7(2) 
di A |. (и bP + a + ho, (к bP + || os 


This formula represents the energy of two noninteracting systems of oscillators 
vibrating with the frequencies w,, and о;у. 


6. The binding energy 


We shall compare now the energies of the ground states for H,,. given by (30) 
and Н,, given by (58). The energy of the zeropoint vibrations for Н с is 


п ho’ _ wh Vo + c2 
S NC LE 


Ex = 
05: 9 


(59) 
and for H oa 


т 
po PRE (60) 


Because of w = Vo? + (02 < wa + ооу it follows that 
Бас (61) 


osc 
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Thus in the case of interaction between the oscillator systems, the zero-point 
energy decreases by the binding energy W originating in the collective motion of the 
system of electrons and ions interacting according to the Coulomb law. 

This energy may be obtained from the relation 

, 
n'h n'h 
co. (m, + оз) — P = "y оу (62) 

The minus sign points to the fact that we һауе to do with attraction. Remembering 
that 025 << оү and assuming п’ ~ N we get 
n'h Nh 


2 Cy EN We ү (63) 


In table I are shown the numerical values for Li, Na, K (Жи =з Ке к= Ж Жез 1) 
of the energy №, arising from the “long-range“ Coulomb interaction responsible 
for the collective vibrations. The energy is compard with the total binding energy. 


W = 


Table I. Binding energies (in kcal/mol) 


PM Calculsted from (63) in Hartree-Fock in Hartree-Fock approxi- 
approximation | mation with correlation 
Li 39 28 14,6 36,2 
Na 26 10 4,1 24,5 
K 22,6 6,5 —0,7 16,5 


The bond considered caused by the collective motion shows some features of 
the metallic bond, viz. the absence of any directional dependence and simultaneous 
acting on a great many atoms (in this respect the metallic bond is strongly reminiscent 
of the van der Waals' forces). 

I wish to thank Professor В. S. Ingarden for suggesting the problem and valuable 
discussions. I also wish to thank Professor M. Günther and Dr A. Rybarski for helpful 


discussions. 


КРАТКОЕ СОДЕРЖАНИЕ 


Галясевич, Обобщение метода добавочных переменных на бинарные системы 


Показано, что в системе состоящей из электронов и ионов (воздействующих электро- 
статически) элементарные возбуждения бывают двоякие: бозонового и фермионного типа. 
Возбуждения бозонового типа можно описать как систему осцилляторов, колеблющихся 
с частотой Мо? + cs?, где o, является частотой лангмуировских колебаний электронов 
(колебаний электронной плазмы ), а оз частотой лангмуировских колебаний ионов (коле- 
баний ионной) плазмы. Эсли пренебречь взаимным влиянием (куломбовским), происходящим 
между ионами и электронами, то элементарные возбуждения бозонового типа избражаються 
двумя системами осцилляторов, колеблющихся с частотами оу и с. Так как V w? + w? 
<®, + @,, To энергия нулевых, колебаний в случае взаимодействия ионов и электронов, 
будет меньше, чем при упущении этого взаимодействия. Разницу энергий обоих этих ну- 
левых колебаний можно считать мерой энергии связи системы происходящей от «дальнего» 
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куломбовского воздействия (ответственного за возбуждения типа бозонового). Оценена 
эта энергия для Li, Ма, К и сравнена с полной энергией связи, вычисленной другими мето- 
дами и измеренной на опытах. Ряды величин будут те же, причем результаты, полученные 
автором оказались ближе данных экспериментальных, чем результаты, вычисленные 
посредством метода Hartree-Focka. 
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LA CHARGE DE L'IMPULSION EN FONCTION DE LA 
SURTENSION DANS UN COMPTEUR С.М. А CATHODE 
METALLIQUE ET EN VERRE 


Pan S. MicHALAK ET А. ZAWADZKI 
Laboratoire de Physique de l'Université de 1642 
(Manuscrit reçu le 5 septembre 1955) 


La valeur de la charge d'une impulsion dans un compteur G. M. d'un type quelconque 
est d’après la théorie de Wilkinson une fonction de la surtension qui s'exprime approxima- 
tivement par une droite qui change son inclinaison en un certain point. Les auteurs ont 
déterminé le coefficient de proportionalité entre la charge et la surtension et ont étudié 
sa dépendance des paramètres géométriques du compteur. Ce coefficient К = АС, où 
C est la capacité du compteur et À une constante qui dépend uniquement du diamétre 
de l'anode. La valeur de ce coefficieni est la méme pour les compteurs ordinaires (à cathode 
métallique) et pour les compteurs de type Maze (à cathode en verre). La connaissance 
de la valeur de 4 est utile dans l'étude pratique d'un compteur G. M. et surtout d'un 
compteur à cathode en verre. 


Introduction 


Dans la description du compteur à cathode externe Maze (Maze 1946) déduit 
une formule d'ailleurs approximative! pour la valeur du potentiel de la surface interne 
de verre : 
ии РК КА 1 

UT ORE 427, Q) 
où V est le potentiel de l'anode, Ур le potentiel du seuil du palier, N la fréquence 
moyenne des impulsions et K un coefficient possédant la dimension de capacité. 

Or ni Maze ni les auteurs qui s'occupaient plus tard de compteur à cathode externe 
et de la formule de Maze n'ont pas déterminé la valeur exacte de ce coefficient (Bassi 


1 [a forme exacte (Blanc 1954) est la suivante: 


2 
e с К (2) +- | 


ет 7 Vo 


(1 + КЕМ)? Ж (1 + КАМ)“ 
la valeur (V — И,)/И, (1 + КЕМ)" est реше par rapport à l'unité et pratiquement négligeable. 
j (65) 
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et Beretta 1949, Favre et Haenny 1953). C'est seulement dans le travail de Blanc et 
Zyngier (1954) que fut déterminée sa valeur dans un cas particulier. 

Les difficultés liées avec la détermination du coefficient K d'un compteur à cathode 
externe résultent du fait qu'elle exige la connaissance de U et de R. En connaissant 
la dépendance R = f(t) (la valeur de R en fonction de la température) pour le verre 
duquel le compteur est fait, on peut calculer К dans chaque cas particulier. La détermi- 
nation directe de U exige soit l'emploi d'un électrométre joint à la face intérieure de 
la cathode, soit un galvanométre intercalé en série dans le circuit du compteur. Les 
mesures les plus soigneuses ne donnaient pas de résultats reproductibles à cause d'une 
grande sensibilité de la fonction К = f(t) et de la variabilité de К le long du palier; 
une autre difficulté provenait du changement de K en passant d'un compteur à l'autre, 
ce qui démontre que la valeur de K dépend de la géometrie du compteur. 

La connaissance de la valeur de K permet premièrement de connaître la valeur 
de U sans utilisation d'un électrométre ou d'un galvanométre grâce à quoi elle rend 
possible une détermination simple de la longueur réelle du palier du compteur Maze 
ainsi que la part de son allongement par l’,autorégulation“ (Michalak, Mowczan, 
Zawadzki 1954). Cette circonstance est spécialement importante pendant les mesures 
d'une source radioactive (l'efficacité décroit avec l’accroissement de N), parce qu'on 
connait la caractéristique du compteur dans chaque régime de son travail. Deuxiéme- 
ment, la valeur de K nous donne la connaissance de la charge Q de l'impulsion d'un 
compteur. Troisiémement, les mesures exactes de K expriment la dépendance de 
K de la géometrie du compteur et du survoltage. 


Signification de K 


K est approximativement le coefficient de proportionalité entre la valeur de la 
charge de l'impulsion et du survoltage 


Qo Ue) (2) 
ce qui donne pour un compteur à cathode externe 
ОЕ UT eU) (Са 


La théorie de Wilkinson (Wilkinson 1948, 1950, 1952) éxamine la dépendance 
de la charge О des différents paramètres paraissants dans l'effet de la décharge 


V l 
НЕ су) > 
Po 2 1n 2 (3) 
où 1 est la longueur du compteur, г, — rayon de la cathode, г, — rayon de l'anode. 
On voit alors que le terme AW1/2Voln(r,fr,) dans la formule de Wilkinson cor- 
respond au coefficient K de Maze. On peut écrire approximativement? 


* La correction introduite par le terme V/V, n'est pas grande et négligeable dans les mesures 
courantes (sa valeur ne dépasse pas 10% en général). 
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K E NE 


ОП 


Ta 


(4) 


où C est la capacité du compteur. Il paraît que la détermination de K pour un compteur 
à cathode interne ne doit pas être difficile, on reçoit pourtant une valeur qui n’est 
pas constante le long du palier. Cette propriété de К est liée avec l'effet décrit par 
Wilkinson, notamment l'amplitude de l'impulsion suivant sa théorie n'est pas pro- 
portionnelle à la surtension tout le long du palier. La fonction О = ДИ — Vo) s'exprime 
par une droite qui passe par le point anguleux m(Q,,. V,,) en diminuant son inclinaison 
approximativement deux fois. Selon Wilkinson la charge contenue dans la gaine des 
ions positifs devient égale au point m à la charge électrostatique portée par l'anode: 
Qn = Ош), Cette théorie explique alors pourquoi on ne peut pas trouver une 
valeur unique de K qui correspondrait à tout le palier. Pour obtenir des résultats 
dépourvus des ambiguités il faut distinguer deux valeurs de K définies de la facon 
suivante 


й: ne Vg) POUT Е И, 
0 

(5) 

2% Q= pU - ине E OV pr 


où V, est le potentiel de la fin РЕ palier, Ку et К, sont liés par la relation 


1 
K, = 9 NES 
Mesures 


Pour effectuer les mesures nous avons utilisé le montage représenté par la fig. 1. 
Les condensateurs „Philips“ avaient une résistance d'environ 1080 pour chaque batterie 
de 160 uF. L’élément RC indique schématiquement la résistance et la capacité du verre 
dans les cas où on a effectué des mesures avec les compteurs à cathode externe. Le 
galvanomètre était du type portable possédant une sensibilité de 6,5 nA sur 1 mm 
de déviation du spot. L’électrométre jouait un réle secondaire: dans une partie seule- 
ment des mesures il servait pour vérifier les valeurs de U (ce qu’exigeait l’utilisation 
du compteur argenté à l’intérieur). Nous avons étudié de différents types des comp- 
teurs à cathode interne et externe dont les longueurs variaient de 11 cm à 51 cm, 
les diamètres cathodiques variaient de 2,5 cm à 5,15 cm, et qui possédaient des anodes 
de 0,05; 0,1; 0,2; 0,25 et 0,3 mm de diamètre. Tous les compteurs étaient remplis 
d’un mélange de 110 mm Hg d’argon et de 15 mm Hg du -méthylal. On a effectué en 


tout 85 mesures. j 


е 


Га ргетіёге série des mesures se rapportait aux compteurs à cathode interne. 
En с dans chaque cas Ја capacité du compteur оп а obtenu pour des tensions 
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pa E ee 


V < V, les valeurs de A (dispersées normalement dans un intervalle des erreurs . 
statistiques) constante pour chaque fil anodique de diamétre donné. D'autre part on 
observait une croissance de А avec l'épaisseur du fil nommé. Il s'ensuit qu'on peut 
fixer une valeur de 4 pour les compteurs possédant une anode d'épaisseur donnée 
indépendamment des longueurs et des diamètres des tubes. 

Nous donnons les résultats numériques dans la suite ensemble avec les données 
concernant les eompteurs à cathode externe. 


* 500 - 5000 V 


Fig. 1. Montage employé à l'étude des compteurs G. M. L'élément RC marque schématiquement la 
résistance et la capacité du verre d'un compteur Maze. 


Les mesures avec les compteurs à cathode externe sont plus délicates. Il faut 
déterminer la fonction R — /(t), puis soutenir thermostatiquement une température 
constante du compteur pendant la durée de mesure, enfin augmenter l'intervalle des 
tensions utilisées (l'effet d’,,autorégulation’’). On a effectué des mesures dans l'intervalle 
des températures entre 17°C et 50°C en utilisant dans ce but des tubesconstruites avec 
des verres de différentes sortes et dont la résistance variait de 140 МО à 2630 mQ. 

Les résultats étaient inscrits dans un tableau de la forme suivante 


ы ve КР IR y-U мее Гота гух 


Г est l'intensité du courant , longitudinal” qui traverse le verre, suivant les génératrices 
du cylindre, de l’anode vers la cathode. Cette correction joue son rôle quand la valeur 
de Г est petite, et on peut la négliger dans le calcul de U = R I. On mesure Г en 
utilisant un compteur rempli d’air sous pression ordinaire et qui par suite ne fonctionne 
pas. La valeur correspondante de I est alors zéro. 

Les résultats obtenus avec les compteurs du type Maze correspondent bien à ceux 
qui sont obtenus avec les compteurs à cathode métallique. La concordance des résultats 


prouve que la théorie de Wilkinson s'applique aussi | bien aux compteurs à cathode 
externe qu'aux compteurs ordinaires. 


tien dni did 
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Les mesures ont été effectuées sous diverses irradiations d'une source radioactive 
(N allait de quelques impulsions par seconde à quelques centaines). 

La figure 2 représente les valeurs de la charge de l'impulsion d'un compteur en 
fonction du potentiel V de l'anode (les cercles blancs, correspondant à l'échelle infé- 
rieure) et du potentiel V — U corrigé par rapport au potentiel interne U de la cathode 


(les cercles noirs, l'échelle supérieure). Le potentiel И est évidemment le méme dans 
les deux cas. 


50 1100 1200 1300 1400 1500 
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Fig. 2. Charge de l'impulsion en fonction du potentiel V (cercles blancs) et du potentiel corrigé V — U 
(cercles noirs). | 
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Fig. 3. La fonction О = f(V — U) Fig. 4. Valeurs du coefficient А pour divers 
ч l diamètres 2 г, de l’anode. 


La figure 3 représente la fonction Q = ДИ — U) d'un autre compteur avec une 
partie plus prolongée du coté V — U dépassant V,, On voit bien que les шш де crois- 
sance de la charge correspondent aux prévisions de la théorie de Wilkinson, 2 inclinaison 
de la seconde partie de la droite expérimentale par rapport à la première est 0,47. 
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МЕЗА ОСОБЕ ВЕ ae E AT AES 


Dans ce cas У, = 1280 volts et Q,, est voisin de Qc, ,,, la valeur de Q,, étant 
2,3 -10-?C et la capacité 1,8 pF. 
Les mesures les plus exactes ont été faites pour l'anode de diamétre 0,2 mm. 


La valeur obtenue pour À est 
А = 7,3 + 0,1 


en exprimant la charge en nC (107°C) et la capacité en pF (10 ??F). 
La figure 4 représente le changement de A avec l'épaisseur de l'anode 


2 r, (mm) 0,05 0,1 0,2 0,25 0,3 
А 6,4 + 0,2 6,7 + 0,2 73401 8,8 + 0,3 9,9 + 0,9 


Applications 


La concordance des résultats obtenus pour différents types des compteurs ayant 
l'anode du même diamètre permet d'utiliser la valeur de À à divers buts pratiques. 
Par exemple la conaissance de À rend possible l'étude d'un compteur Maze sans 
savoir quelle est la résistance du verre et en n'utilisant que l'oscillographe cathodique 
(Zawadzki 1954) (on n'a pas besoin d'échelle ni de galvanomètre). Comme au point m: 
О„ = Qu ," nous pouvons écrire deux équations, contenant deux inconnues O) 
еи: 
Qo» = Vin — Un) (6) 
V тъ Un 


0 


Qn= АС Vin — Un = о). ^ | (7) 
L'étude des impulsions à l'aide d'un oscillographe nous donne la valeur Z,, (montrée 
par le voltmétre de la source de tension) et aussi la valeur V, (par prolongement de 
la droite jusqu'à son intersection avec l'axe des abscisses). 

La figure 5 représente la dépendance entre l'amplitude des impulsions (dont 
la lecture était faite sur l'écran d'un oscillographe au moyen de papier millimetré 
iransparent et d'une loupe) et la tension appliquée au compteur. 

Les mesures fournies par l'oscillographe sont évidemment relatives à cause de 
ne pas pouvoir passer directement de la hauteur de l'amplitude lue à l'oscillographe 
à la valeur réelle de l'impulsion produite par le compteur (effet des capacités de liaison 
et de la constante de temps du couplage). Il s'ensuit à fortiori que l'on ne peut rien 
savoir sur la valeur de la charge créée dans le compteur. On voit d'ailleurs que la con- 
naissance de cette valeur n'est pas indispensable pour faire nos calculs. On peut pour- 


tant trouver d'une facon indirecte la valeur de la charge produite, à partir des relations 
. (6) (7) on obtient 


V, i ык А 


yar 
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En connaissant la valeur de (7, et les paramètres géometriques du compteur on obtient 
facilement Q,,. D'autre part en connaissant la valeur de U,, on peut déterminer sur 
un graphique corrigé (les données obtenues directement expriment l’amplitude en 
fonction de V et il nous faut passer à V — Unm) la valeur de О pour une surtension 
quelconque. Si on dispose d'un scale, on peut déterminer aussi la résistance du verre: 


0, 
Om N —I— тъ. (9) 


Dans le cas d'un compteur métallique on n'a méme pas besoin de prolonger les mesures 
jusqu'au point de disconinuité et il suffit de trouver Бу. Il va sans dire que l'étude 
сотріёе d'un compteur C. M. ordinaire peut s'effectuer en n'utilisant qu'un oscillo- 
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Fig. 5. Amplitudes des impulsions en fonction du potentiel de l’anode (non corrigé par rapport à U) 
lues sur l'écran d'un oscillographe. L'échelle de l'axe des ordonnées est arbitraire. 


graphe et une source de tension munie d'un voltmétre, ce qui facilite appréciablement 
le contróle des compteurs dans le laboratoire. La courbe rebroussée s'applatit visible- 

: : е j mor 
ment après avoir dépassé la valeur V, (la fin du palier) ce qui permet d'établir directe- 
ment quelle est la longueur du palier. 


Comparaison des résultats obtenus avec ceux de Blanc et Zyngier 


Blanc et Zyngier (1954) ont trouvé K = 12 pour un compteur d'une longueur 
| = 16,4 cm de 2r, = 0,1 mm et de г, = 1 ст. La capacité de ce compteur est alors 
C = 1,72 pF ce qui donne A = 12/1,72 = 7, la valeur de À trouvée par nous pour 
une anode de méme diamètre est 6,7 + 0,2. Leur compteur était rempli De mélange 
de 94,5 mm Hg d'argon et 10,5 mm Hg d'alcool (une proportion pro 50% 
plus petite de la vapeur organique et une pression totale de 20 mm Hg inférieure que 
chez nous) et malgré ça il y a une concordance des résultats. du peut alors supposer 
que la valeur de À ne dépend probablement pas d'une facon appreciable de la composi- 
tion et de la pression du mélange gazeux au moins dans les limites ordinairement 


adoptées. 
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КРАТКОЕ СОДЕРЖАНИЕ 


Михаляк и Завадский. Заряд импульса в функции перенапряжения в счётчике С. М. 
с металлическим и стеклянным катодом. 


Величина заряда импульса в счётчике С. М. любого типа является, согласно теории 
Вилькинсона, функцией перенапряжения, приблизительно выражающейся посредством 
прямой, меняющей свое наклонение в известной точке. Авторы определили коэффициент 
пропорционального взаимоотношения между зарядом и перенапряжением и изучили его 
зависимость от геометрических параметров счётчика. Этот коэффициент К = AC, причём 
С является ёмкостью счётчика и А постоянной, зависящей исключительно от диаметра 
анода. Значение это коэффициента такое же для обыкновенных счётчиков (с металличес- 
ким катодом), как и для счётчиков типа «Маза» (с катодом стеклянным). Знание значения 
À полезно при практическом изучении счётчика С. M., а в особенности счётчика со стек- 
лянным катодом. 
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NOTE ON THE SLI (д,«) “НЕ REACTION 


Ву JERZY SAWICKI 
Institute of Theoretical Physics, Warsaw 
(Received July 7, 1955) 


Resnick and Inglis (1949) have given the compound nucleus theory for the 
SLi (d,a) “Не reaction and have compared their results with existing experiments up 
to energy 3.75 Mev. In this energy range satisfactory agreement with experiment was 
obtained for the angular distribution. | 

However, for much higher energies the pick-up process of the odd proton and 
neutron loosely bound with the «-core in Ше Li nucleus by the incident deuteron seems 
to be most important. 

Recently (Dabrowki and Sawicki 1955) a theory was given for the 9Li(n,t) “Не 
reaction based on the simple ,,a-particle-deuteron model“ of the 6Li nucleus. Reason- 
able agreement with experiment was obtained. For the (d,a) reaction similarly as 
for the (п,)) reaction the pick-up process is most probable when the odd n-p system 
in the Li nucleus is bound. The less important Coulomb effects are neglected. 

Taking into account the indistinguishableness of the a-particles we may obtain 
the Born approximation differential cross-section in Ше c.m. system similarly as it 
has been done for the (n, t) reaction (Dabrowki and Sawicki 1955) in the form; 


do Ma Ме ke 1 KO + Кл —0)? а 

dQ (ол!) js Dy Me Е i: 2, (1) 
where М = (3/2) M and М *, = 2M are respectively the reduced masses of the 
incident deuterop and tho *pick-up^ a-particle, ka and ky are their respective 
momenta in the c. m. system, ии and ив are the incident deuteron and the ‘Li 


nucleus magnetic numbers, 0 = — (kg, Ка), 
J 10) du Le (ДЕТ Q) 
where Xa X, (Е) is the internal wave function of the a-particle in the Li nucleus, 


Ха E Solo, 


n? 
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pi On др) Har (ts 0, м) = Sol Ou 095 On Ор) Ra (ust) — the antisymme- 


(71) 
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trized function of Ше, *pick-up^ a-particle, *) 
та т Tm — 
ge DE и ende Р , w = اگ‎ 
2 2 y2 y2 
are the internal coordinates of the *pick-up^ a-particle, 7 is the vector from the «- 
particle to the center of mass of Ше n-p system in the Li nucleus, V = V(r) is the 


“q—d* interaction in the ‘Li nucleus, y, is defined by 
Yo = Xa (Est. V23.,,0,) Sra, (ало Op) р (УЗие Fa, (3) 


where the o's are the respective spin variables of the four nucleons, x, is the internal 
wave function of the SLi nucleus, Sug, (др б,,) is the incident deuteron spin function 


Га Us bos Tns — Ре 


and y(V2w) its internal wave function. According to our “a + d“. model х; may 
be approximated in the form 


Xs = Ха (5)0)5,, (9,05) v( V2 v). (4) 
where Ф (г) is the S—state “a + 4“ system internal function. 
Consequently we have 


1 1 
= У |0) + aol = =5 |, OL (8) LG-91G-95 — 6 


where À 
I, = f dudvdw е R,(u,v,w) y (V2) v (V2 ш) (5a) 
I, = Í dre Ир) Ф() (5b) 
ме 
and а = e Кы ky p= СА — Ка, - 


Probably the angular distribution is little sensitive to the choice of the function 
Ra similarly as for the (4,6) reactions to the choice of the triton function R, (Newns 


do 
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1952). For the sake of the calculations the gaussian function given by Bruno (1948) 
was used. For this function the angular distribution does not depend on the choice 
of the deuteron function (о). 

For V(r) we use a square well with Ry = 6x 10-13 cm — the better for higher 
energies (Dabrowski and Sawicki 1955). The angular distribution of the «-particles 
(symmetrical about 6 = 90°) is given in Fig. I for the incident deuteron energy 20 Mev. 
The second big maximum at 6~40° is due to the great О value (24.49 Mev) and the 
square well «—d interaction. Unfortunately no experimental data are at present 


available. 


The absolute values of E: for the energy 20 Mev are respectively 0.2 
9-0 

X 1073? cm? for y(o) being Ше usual square well function and 0.4 х 10-39 cm? for y(o)- 

hulthen (e.g. used by Kruse, Malenka and Ramsey (1953)). The reduction of these 


values due to the Coulomb effects is estimated as small. 
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COAL SURFACE ANALYSIS BY MEANS OF AN ELECTRON MICROSCOPE 


Ву ANTONI FELTYNOWSKI AND IRENA GLASS 
Institute of Physics of the Polish Academy of Science, Warsaw 


(Received September 26, 1955) 


It is the purpose of this paper to elaborate а method for the examination of 
pitcoal surface by means of an electron microscope. Coal surface has been examined 
in an electron microscope already in 1943 by Preston and Cuckow, yet even now 
the obtaining of exact replicas is still a difficult problem. In our investigations we used 
chromium shadowed aluminium replicas and it seems that we have obtained better 
results than previous authours, which have been using negative formvar replicas and 
positive polystyrene replicas. 

Samples of coal in the form of cubes with edges of 2 cm have been mechanically 
polished by the method of Stach (1928) and Benes and Růžička (1952). The surfaces 
thus prepared have been submitted to the action of 10% chromic acid during 10— 
20 seconds, washed with warm water and after being dried covered with saponlack 
(nitrocellulose 10%, tricresil phosphate 2%, butyl acetate 88%). Replicas of samples 
not etched but only washed with ethyl alcohol were also made. After 36 hours the 
saponlack layer was detached by means of a paper band covered with a layer of glue. 
At first the saponlack layer was incised with a scalpel at a distance of about 2 mm from 
the edge of the cube. The paper band was stuck to the edge (with the incision) 
of the saponlack layer. By means of the band the saponlack layer was separated 
from the coal surface in the form of the pieces free from the paper. The pieces 
of the saponlack replica were attached to a miroscope glas-slide outside the side 
which before touched the coal. The replica on the glass-slide was first shadowed 
with chromium in vacuo at an angle of 15°, and afterwards it was covered with 
a thin layer of aluminium (of about 50 À) by direct evaporation. The replica, 
after being detached from the glass-slide with a scalpel, was placed on a phosphor 
bronze grid, which usually serves to support the specimen in the electron micro- 
scope. It was placed in such a way that the saponlack layer was at the top and the 
aluminium layer touched the grid. The grid with the replica was moistened with amyl 
acetate and afterwards immersed in it for 12 hours. The saponlack layer dissolves 


(15) 
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then completely and the aluminium replica shadowed with chromium either remains 
on the grid or rises to the surface of the liquid. In the second case it should be 
carefully taken out with the grid. After being dried the specimen is ready for exa- 
mination in the electron microscope. 

By means of this method we have begun investigations of the microstructure of | 
coal strata for the Mining Institute. The inserted риол illustrate the results 
obtained. 
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Fig. 1. Zabrze-East coal, unetched, magnification 5000. 
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Fig. 3. Zabrze-East coal, etched with chromic acid, magnification 4500. 
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dig 4 Zab ERR coal, etched with chromic acid, magnification 4500. 


- Fig. 5. Kennel coal, magnification 4800. 
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| | ; Fig. 6. Kennel coal, magnification 5000. 
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